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NOMENCLATURE 


A 

a 

a , b » ^ 



c 


c* 

c 

P 

D 

D 

£ 

m 

F 




area 

local sound speed 

empirical spray coefficients, Eq. (2) 
drag coefficient 
thrust coefficient 
approximate evaporation coeff: nt 

nozzle discharge coefficient 
mixture ratio 

characteristic exhaust velocity 

specific heat at constant pressure 

droplet diameter 

mass median droplet diameter 

Rupe mixing efficiency factor (Ref. 7) 

drag force 

gravitational coefficient 

specific impulse 

heat of vaporization 

thermal conductivity 

droplet evaporation coefficient 

molecular weight 

droplet mass or spray size group mass 

rate of change of mass 

droplet concentration (no/volume) 

number flowrate of droplets 

number of injection elements 

Nusselt number 

number of stream tubes 

pressure 


Iv 



u (OH 


p* 

Be 


V 

"ooi 

w 

g 

X 


x»y.2 


Z 


Z 


0 


= pressure corresponding to sonic flow 

» Reynolds number 

» universal gas constant 

» nozzle throat radius ratio (curvature/opening) 

* nozzle throat opening radius 

= radial coordinate 

* stream tube path length 

- temperature 

= velocity 

* propellant mass flux at a spatial mesh point 

= propellant mass flux contribution from an injection 

element to a mesh point 

“ flowrate 

» w(x,y,z) for x» 3 r» 0 , z»l 

* stream tube gas flowrate 

* non dimensional distance from nozzle throat 

* rectangular coordinates (referenced to an injection 

element origin) 

« axial coordinate 

= initial plane for beginning spray combustion analysis 


GREEK LETTERS 


a 

a,6.V 




c 

c 

e 

e 

"At 

n 

0 

u 

p 

CT 


= nozzle wall angle (from chamber axis) 

» element relational angles (see page 72) 

» ratio of specific heats or adiabatic expansion coefficient 

• pseudo-impingement points for gaseous and liquid propellants 

for a gas/liquid injection element 

» chamber contraction ratio 

• nozzle expansion ratio 

■ decimal tolerance, convergence on throat area 
“ efficiency factor 

“ angular coordinate 

• viscosity 

• density 

■ surface tension 
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SUPERSCRIPTS 


( ) 
n 

prime, ( )' 


average value or concerned with one-dimensional solution 

t^i 

concerned with the n droplet size group 
value reduced by evaporation 


SUBSCRIPTS 

c 

cr 

d 

E 

e 

f 

g 

i 

inj 

j 

i 

o 

0,0 

s 

STC 

T 

TDK 

v,vap 


chamber 

critial point property 

droplet 

element 

expansion section 
fuel, spray fan 
combustion gas 

th 

concerned with i stream tube, summation index 
inj ection 

concerned with propellant 

liquid 

oxidizer 

initial or stagnation value 
surface or stream tube path 
stream tube combustion model 
throat 

two-dimensional kinetic model 
vapor or vaporization 
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1. INTRODUCTION 
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The Distributed Energy Release (DER) combustion analysis computer program 
described this report has been developed to fill a need in, and to become 
a part of, the JANNAF system of related computer programs for performing bi- 
propellant liquid rocket performance analyses. The DER program is concerned 
with formation, distribution, flow and combustion of liquid sprays and com- 
bustion product gases in conventional (cylindrical) rocket combustion chambers. 
The main spray flow direction must be essentially parallel to the chamber axis, 
but the injector surface can deviate from a simple flat plate normal to the 
chamber axis. 

Although an earlier evaporation coefficient version of DER (Ref.l and 2) also 
had a capability to analyze non- equilibrium, two-dimensional supersonic exhaust 
nozzle flow, this December 1971 version simply provides output data for subse- 
quent input to another computer program (TDK) for performing that analysis. 
Major modifications included in the current version are: 

1. Simplified injection element specification for inputting the LISP program 
block, 

2. Replacement of the spray flux and mixture ratio distribution plotting sub- 
programs with simpler, reduced-run- time routines, 

3. Provision of a partial LISP capability for analyzing bipropellant injection 
involving one gaseous propellant, 

4. Removal of a modified, long-form- option revision of the TDK program block 
from DER, 

5. Provision of a subroutine for punching out nozzle throat-plane data in the 
form of a partial data deck for subsequent use in running cm improved 
multi-stream tube TDK program. 

6. Interface option to punch cards In LISP for 3DC and to read punched 
cards from 3DC in STC. 

Regarding engine performance calculation, DER computes only characteristic 
velocity efficiency; all other performance parameters are presumed to result 
from the subsequent TDK program analysis. 

A technical description of the DER computer program is given in Section 2 of 
this report. Beginning with Section 3, the remainder of this report is 
concerned with DER computer program usage. Section 3.1 shows the recommended 
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program overlay structure. Section 3.2 contains a brief statement about each 
of DER's subroutines, indicating the major functions that each performs. 

Section 3.3 is concerned with the provision of the rather considerable amount 
of input data needed to make DER computer program runs. Sample data code sheets 
show the correct sequence and formats for variable values and include notes 
concerning options to bypass the input of certain data. A separate list defines 
the variables and indicates their required dimensions. Section 3.4 contains 
some notes concerning program operation, e.g., core size requirements, output 
limits and execution times. Section 3.5 describes the conputer program output. 
Literature cited in the body of the report is listed in Section 4 and correla- 
tion coefficient values used in the LISP subprogram are tabulated in Ref. 6. 
Finally, appended to the report, as a separate volume, are a listing of the 
computer program (Rev. Jan. 74), an exanple data deck listing and portions of 
typical output data. 


la 
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2. COMPUTER MODEL ANALYSIS 

Basically, two major subprogram blocks comprise the DER computer proff'' ' 

These subprogram blocks can be executed Individually or fully coup. as a •r;le 
program. The first is the Liquid Injection Spray Pattern (LISP) projiram. Li. "^ 
utilizes correlations for the propellant mass flux distributions for single- 
injection elements and, based on element design, location and orientation data, 
combines outputs from all injection elements to calculate a complete mass flux 
distribution in a plane normal to the combustor's axis and a distance z^ down- 
stream of the injector. Other outputs from this program block include mean pro- 
pellant droplet sizes and velocity vectors, approximate calculations of the per- 
centages of propellants gasified upstream of z^, and data regarding the resultant 
combustion gas flowing through the plane z^. Cylindrical combustor coordinates 
are used, 'fhe output from LISP provides the necessary description of the two- 
phase flowfield for initializing the second major subprogram block, the Stream 
Tube Combustion (STC) program, STC is designed to accept propellant flux data 
at r, 0, Zq mesh points from LISP, combine them into several axisymmetric stream 
tubes and analyze spray gasification and combustion in all stream tubes 
simultaneously as the flow progresses downstream of Zq, STC's basic method 
involves finite difference solution of coupled ordinary differential equations 
by marching along the chamber length. 

An optional procedure Is provided to supplement the DET analvsls by executing 

i 

a 2-phase, 3“dImensIonal combustion program, 3DC, In between the LISP and STC 
subprogram analyses. In the normal DER analysis, spray Is considered to travel 
In straight flow paths completely independent of gas flow paths throughout the 
region modeled by the LISP subprogram and then instantly forced at the start of 
the STC analysis to change direction to paths coincident with gas streamline 
flow. The 3DC analysis physically models the mixing of spray and gas as the 
original spray flow paths are redirected by the transverse drag force with the 
gas flow. 

*This is the 3C-C0MBUST (3DC model) reported on in Ref. 5, only with slight 
modifications in order to accept punched card output from LISP and to provide 
punched card output for STC. 
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At some distance upstream of the combustor exhaust nozzle ’r throat, the STC 
program block calls upon subprogram block TRANS to provide information about 
the transonic flow's spatial pressure distribution. The geometric relation- 
ships among LISP, STC and TRANS program block analyses are indicated in Fig. 4 
(page 23). Application of the Two-Dimensional Kinetic (TDK) nozzle analysis 
computer program (Ref. 3) is also indicated there, although this is no longer 
an integral part of DER but must be accomplished by subsequent use of DER out- 
put data. It is because TRANS and TDK are formulated in two dimensions (r, z) 
that STC is designed to assemble asixymmetric stream tubes at plane z^. 

Descriptions of the analytical bases for the major subprogram blocks are given 
in this section. 



2.1 INJECTION AND ATOMIZATION: 

LISP SUBPROGRAM BLOCK 

The LISP computer program block has been developed , applied and improved in 
several successive stages. It vas initially structured as the first of several 
individual programs in a system of computer programs for analyzing injector/ 
chamber compatibility lor liquid/liquid bipropellant combinations (Ref. 4). 

That version was adapted to form a LISP subprogram block in the initial DER 
computer program (Ref. l). Subs j^uently, the LISP input data requirements 
were simplified, subroutines for producing computer-plotted mass flux contour 
and profile graphs were simplified and the formulation was expanded to include 
coaxial jet injection of gaseous fuel/liquid oxidizer propellants (Ref. 5). 
Similarly, gas/liquid triplet injection was added during the current DER program 
effort. 

The LISP computer program calculates spray mass fluxes at mesh points (r,0,Zg) 
by straightforward summation of mar ^ fluxes from all individual injector 
elements 

^ ^EL ^ 

V(r,6,z ) * 2 w (r,e,z ) (l) 

'^ 1=1 '' 

The method can be o ed if: (l) the individual injector elements have predic- 

table spray flux patterns which can be measured €uid correlated with the main 
direction of spray flow parallel to the chamber axis; (2) the individual 
spray patterns of the various elements are not destroyed between the element 
impingement points and the pleme by collisions between neighboring fans; 
and (3) vaporization of injected spray mass between the element impingement 
and plane z^ can be estimated. A detailed accounting is not made for any back- 
flow of spray toward the injector; overall axial flow continuity is enforced 
at the LISP collection plane* 

The necessary correlation for w^ (r.QfZ^) for use in Eq. (l) is based upon the 
shape of single element spray flux distributions determined from cold-flow 
spray measurements. Flux patterns were fitted to the generalized expression: 
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v.(x,y,z) - ^ I [l + Cj(^) + C 2 (^)^ + [c^C*) 


♦ C,(i)^fl >C (^) I- 

4'z' j L 5 * O * J J 


-a(!) -be) 


( 2 ) 


which is applied separately to each propellant frora an element. The (x,y, 2 ) 
coodinate system in Eq. (2) is referenced to the element's impingement point 
from which its spray is presumed to emanate. The empirical coefficients 
a,h,w^, and C through are functions of such pai,uneters as element type 
(doublet, triplet, etc.) impingement angle, orifice diameter, impinging stream 
momenta, orifice length, and manifold effects. The coefficients are obtained by 
correlating data from injector cold-flow simulation tests, in which immiscible 
propellant simulants are collected in a plane beneath the injector and measured. 
The form of Eq. (2) was chosen because it satisfies continuity, predicts the 
observed inverse square relationship between mass flux and distance from the 
impingement point, and because closed form integrals of Eq. (2) ai. 1 its x and 
y moments over the x. y plane allow straightforward evaluation of the empirical 
coefficients, using experimental cold-flow data. A detailed discussion of the 
coefficient evaluation procedure, including an actual example computation, is 
given in Ref. 6. That document is primarily a catalog of the correlation 
coefficients which are currently in the LISP library of coefficient subroutines. 

LISP transforms the x,y,z coordinate systems of the individual injection 
elements to the r, 6,z coordinates of the combustor and then suras the propellant 
mass fluxes from the individual elements to each r,6 mesh point in plane z^. 
Advantage can be taken of injector design symmetry to define a mesh system which 
occupies only a fraction of the chamber. This reduces the amount of input data 
and the LISP program execution times, but care must be exercised in defining 
the injection elements, to ensure that the mesh system receives the correct 
geometric proportion of the total propellant injection rates. 

LISP assumes that droplets travel as rays from the element orir >ns (e.g., 
impingement points) to the mesh points and that injection kinetic energy is 
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conserved. The velocity of a given drop is therefore given by its injection 
velocity and the direction cosine of its path from tlie injector impingement 
point to the mesh point. When expressed in a common reference frame (the 
chamber geometry) the mass-veighted drop velocity vectors from individual 
elements can be summed at a mesh point to give mean velocity vectors ”dr* 


and u 


de* 


Propellant fluxes from neighboring elements are considered to pass through one 
another without interference. For many rocket designs, with moderately small 
elements or moderately widely-spaced elements, this assumption appears to be 
reasonably valid. Wlien a spray was examined from a non-continuum, kinetic 

theory viewpoint, the mean-free-path of a drop was estimated to be about l/2 to 

. 2 

1-inch for a typical rocket injection condition of about 1 Ib/sec-in . 

Operation at higher injection densities (higher chamber pressures and/or lower 
contraction ratios) would tend to invalidate this assumption, but no limit con- 
ditions have been established. 


Gaseous propellant injection is currently treated in LISP as if it were a low 
density, non-continuum spray. Modifications to the LISP logic were thus 
minimized, but this treatsMnt obviously limits the program to analysing aiagle* 
elsmanit injectors or mnltlple-eleatent injectors whose elements are tar enough 

apart that their gas flows do not interact strongly. With cloaely-spaced 
elements, the gas distribution can be rather poorly predicted; a symptom of 

this condition may be the calculation of axial gas velocities between elements 
that are comparable with those downstream of the elements. 

A very essential part of the combustion field initialization performed by 
LISP is the assignment of propellant droplet size distributions. In the DER 
computer program, LISP computations are concerned only with a mass median 
diameter (D) for each propellant's spray. Later, during STC program block’s 
initialization of stream tubes, the sprays are distributed into a discrete 
number of droplet size groups. The magnitudes of the D’s frequently have a 
direct, strong influence on the propellant combustion efficiency computed by 
DER. Thus, an effort should be made to ensure that the most realistic values of 
D are supplied to LISP. 

For some injector types, mean droplet sizes are calculated within LISP if they 
are not supplied by the program user. Here, empirical equations are used which 
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relate D to injection element hole sizes and injection velocities. The basic 
formulae were derived from cold-flow studies of molten wo.« jets; empirical 
adjustments are included to account approximately for differences in liquid 
and gas properties between those experimental and combustor conditions. A 
discussion of LlSP's D formulae and guidelines for selecting D's to input 
to LISP are given in Ref. 6. 


Partial propellant evaporatio.n upstream of is calculated by a simplified, 
integrated evaporation expression 


w'(r,e,ZQ) * w(r,e,ZQ) 


1 


'ill 

REPRODUCniTT.rrV^ OF THE O) 
ORIGINAL PAGE IS POOR 


where w* is the liquid spray flux actually arriving at the point (r, 9, z^). 

The coefficient is related to the evaporation coefficient k’ used in the 
subsequent spray combustion analysis* However, beca’ise the liquid sprays are 
not fully atomized over the entire Az distance, values of , including a con- 

vections Nusselt number, are usually assumed to be only about l/5 to l/4 of the 
stagnant values of k* • The propellant vapors said to be generated by this 
calculation are summed over all mesh points to yield a single overall vapor flow 
rate for each propellant* Use of such a simplified evaporation expression is, 
to some extent, justified by the relatively small percentage of evaporation in 
the spray formation zone. 


To ensure continuity of both propellants, the LISP program performs summations 
over all mesh points of local propellant flow rates and then scales each pro- 
pellant's mesh point flows and fluxes by the ratio of its total injected flow 
to its mesh point sum. This normalization process thus corrects the total 
flows for computational simplifications ani errors. 


Input to the LISP computer program consists of the number, location orientation, 
size, geometry, and type of injector elements, together with the geometry of 
combustion zone mesh network and general data concerning the propellant densities 
and injector pressure drops. Up to 60 injector elements can be considered* As 
many as 40^3 combustion zone mesh points can be prescribed* 
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2.1A 2-PHASE, 5-DIMENSIONAL SUBSONIC FLOW: 3DC INTERFACE WITH DER 

COMPUTER PROGRAM 

Two distinct models for analyzing flow in the combustion chamber exist in the 
DER coiaputer programs: one is the LISP model which analyzes a region starting 

at the injector face and the other is the STC model which analyzes the region 
downstream of the LISP region. These models are integrated in DER to be 
executed sequentially in a single computer run. However, by executing the two 
DER models separ. ely, the overall analysis can be expanded to include the 3DC 
computer program for modeling 2-phase, 3-dimensional flow in an intermediate 
region. Interfacing of the DER models with 3DC is not fully automated. To 
assist the user in running these models in sequence, the bul’: of the flow data 
which must pass between them is transmitted via computer-generated punched data 
cards: LISP generates data cards for 3DC, 3DC accepts these cards and punches 

cards for STC, and STC reads the 3DC-generated cards. 

Inclusion of 3DC imposes a severe requirement on the mesh size used in LISP. 

The sector which is representative of the cross section must contain seven 
angular mesh positions and 15 or less radial mesh positions. LISP can accomo- 
date a mesh size up to 20 by 20 when 3DC is not used, but must have the same size 
mesh system as 3DC when the latter is used. For som<^ injector designs, the 5DC 
mesh size may be too course to obtain satisfactory ^ t-formance predictions. 

The length of the LISP region is especially critical in detemiining the propel- 
lant mixing efficiency when the 3DC analysis is bypassed, because the length 
controls the amount of overlap between neighboring spray fans. Ideally, this 
length is not critical when it is kept short and the spray spreading is handled 
by the 30C analysis. However, due to practical considerations of modeling 
spray flow with a finite element representation, the LISP region length may or 
may not be very critical, depending on the particular spray-spreading character- 
istics. In the LISP analysis, spray at any mesh point is moving in as many 
directions as there are injection elements. In interfacing with 3DC, a single, 
mean spray velocity vector replaces the multiple spreading spray velocity 
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vectors at each mesh point. Therefore, the modeling of spray spreading using 
30C with LISP loses some validity when propellant Is Injected from more than 
one Injection element. In general, user judgment Is required to select an 
appropriate length for the LISP region, whether 30C Is used or not; however. 
In some cases, the judgmental selection Is not critical. 
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2.2 SUBSONIC FLOW AND COMBUSTION: 
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STC SUBPROGRAM BLOCK 

Selected data computed by LISP are transferred (via scratch data unit) to the 
STC subprogram block. There, by one of three alternate methods, several mesh 
points' propellant fluxes and flow areas are combined to form one of the stream 
tube flows to be analyzed by STC.* Model solutions for spray gasification 
and combustion are obtained numerically for several systems (one for each 
stream tube) of simultaneous ordinary differential and algebraic equations 
by starting from known conditions at the LISP collection plane and marching 
downstream in small axial steps. 

Input to the STC computer program consists of chamber wall profile, propellant 
properties, combustion gas properties and either (1) initial-plane gaseous 
flowrate and mixture ratio and spray flowrates, velocities, and droplet diameters 
for all spray size-groups entering each stream tube or (2) data from LISP from 
which these variables can be calculated. Up to 40 stream tubes can be 
initialized with as many as 12 spray size-groups (fuel and oxidizer combined). 


♦When 3DC is used, these dAta are transferred via punched cards generated 
by 3DC. 
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2»2.1 Stream Tube Initialization from LISP Data 


Data transferred to STC from LISP are: propellant mass fluxes in tlieir Injected 
states, mean droplet velocities and mass median diameters at each mesh point; 
mesh point coordinates; and total initial plane flow and how much of it is 
gasified for each propellant. If there is no spray present for one of the 
propellants, its gaseous mass flux at each mesh point is retained, otherwise, 
simple bulk propellant vapor flows are transferred. In either case, the gas 
mixture ratio is considered to be uniform (constant) across the r,0,z^ plane. 
Axisymmetric stream tube flows may be initialized from these data by means of 
the following options: 


(1) Ail mesh points along each circle (r « constant) of LISP*s mesh points 
are combined into one stream tube. Gasified propellants are handled in one 
of two ways, depending upon whether one propellant is completely gasified or 
not: 

a) If there is no spray of one propellant, the bipropellant gas flux dis- 
tributions to the mesh points computed by LISP may either be retained 
or be averaged out to form a uniform bulk gas flux before the mesh 
points are combined into stream tubes. 

b) If there is spray of both propellants, the gas mass fluxes are 
initially approximated as being uniform: 


V 


"ij 



(4) 


Then the gas mixture ratio at each mesh point is said to be equal to 
the spray mixture ratio there: 


c 


ij 



ij 


(5) 


In general, however, those two assumptions will not be compatible with 
conservation of propellant species flow rates, e.g.: 



4 r 

ij 
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UL 


1 + c 


ij 


( 6 ) 
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Therefore, the fuel and oxidizer contributions to each mesh point's 
gas flow are scaled separately to preserve species continuity: 
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c. . 


ij 




c. .w 
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lil 


1 ^ c 


ij 


(7) 
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These definitions complete the specification of propellant flows at 
each mesh point. 


This initialization method may be appropriate for injectors that form essentially 
axisymmetric flows. When applied to injectors which reduce angular gradients 
in local propellant mixture ratio, however, it can effect substantial mixture 
ratio averaging and result in overcalculation of combustion efficiency. 

(2) Mesh points are combined into stream tubes on the basis of position and 
local mesh point mixture ratio. Distribution of the gaseous propellant flows 
is the same as described in (la) or (lb), above. Following distribution of the 
gases amont the mesh points, a wall boundary layer stream tube is established 
by combining all the mesh points at the wall. If that stream tube does not 
contain more than a specified percentage of the total flow, the next inward 
circle of mesh points will also be combined into it, etc., until it does. 

Then the remaining LISP circles of mesh point are divided into a specified 
few (2 to 4, perhaps) circular or annular zones having roughly ^qual propellant 
flow rates* as illustrated in Fig. la. 


Within each of these zones, the mesh point flows are accumulated into stream 
tubes according to their total propellant mixture ratios, rather than to their 
positions. The number of stream tubes per zone is specified and they are 
assigned roughly equal propellant flow rates. The lowest mixture ratio mesh 
points are combined into the first stream tube until its fraction of the 
zonal flow rate is reached, the next lowest mixture ratio mesh points are 
assigned to the second stream tube, etc., as illustrated in Fig, Ih, 
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a. Gaometric Assignioent of Meah Points to a Wall Boundary 
Layer Flow and to Two Approximate ly-Equal Flow Raie Zones. 



Figure 1. Illustration of Assignment of Mesh Point Flows 
to Stream Tubes of Like Mixture Ratio 
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Finally, the resultant stream tubes are arbitrarily assigned radial positions 
within their respective zones » with the fuel-rich stream tubes lying inside 
of the oxidizer-rich ones. 

This method preserves the angular averaging objected to before o^ily at the 
wall and is accepted for a fraction of the flow in order to get a wall- 
bounding stream tube that is characteristic of the mean wall mixture ratio* 

For the remainder of the flow, the non-physical combining of mesh points on 
the basis of mixture ratio has been found to effect only modest changes in 
calculated mixing efficiencies from those based on the full LISP distributions. 

2.2*2 System of Equations 

The system of equations for the i— stream tube is: 

Gas Phase 
Continuity: 




Momentum: 


J.n 


* F. 


J.n 


Adiabatic Energy Equation: 




01 


where 


■^oi • -y(cj). and - M^(c.) 
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( 11 ) 
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are tabulated* and 
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oi 
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( 12 ) 


This corresponds to frozen expansion to local conditions from stagnation 
equilibrium. Frozen values of the specific heat ratio, y, are used. 


The local stream tube gas mixture ratio is obtained simply by integrating the 
evaporation rates to get gasified flowrates! 


^ “^0 

Mixture Ratio: 


( 13 ) 


( 14 ) 


State: 


0 . 


p M . 

wi 
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Spray Phase droplet size group of propellant) 

Mass Continuity: 

(«a ")- A 1 = - A (m."). 

ds L ^dj 'x dj 'i s^J ' j 'i 


Drop Number Continvity; 

T [ (Nj ")- (uj ")- a 1 

ds L ' dj 'x ' dj 'x sj 


= 0 


( 15 ) 


( 16 ) 


( 17 ) 


These combustion gas properties are obtained from separate calculation of 
equilibrium chamber conditions for several mixture ratios and the nominal 
chamber pressure for a particular case being analyzed. They are also rela- 
tively weak functions of chamber pressure, but this dependence is neglected. 
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or, equivalently 




(Nj.”). (m,.”). a = constant* 
' dj ' dj 'i 8. 


Momentum; 


d_ 

ds 





■ *,[ec 


( 18 ) 




The independent variable in these one-dimensional flow equations is the stream 
tube path length or flow direction, s^. This variable is related to the stream 
tube’s cylindrical (r,z) coordinates through the differential expression 


ds. 

1 



. & 1/2 


( 20 ) 


where r. is the stream tube’s mean radiua. 
1 

solution, however, approximations are used 
wall is parallel to the axis and that 



For numerical stability in the STC 

that ds . “ dz where the chamber 
1 


( 21 ) 


in the nozzle. The basis for Eq. 21 may be seen by examining Fig. 2* 

In this formulation, Ag^ appears as a dependent variable for which a solution 
is to be found. The gas phase equations are constrained, however, in terms 
of z-plane area: 


S (z) = aJz) 

i i 


■*^his equation states that the number flowrate of droplets in each propellant 
spray group size is constant. However, if only a single size group (i.e., a 
monodisperse spray) is specified for c. propellant, the program holds the drop- 
let diameter constant for that size group and its number flowrate is decreased 
as spray vaporization proceeds. 
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Figure 2. Schematic Illustration ol Assor 
Conical Convergence of Stream 1 


Therefore, the foregoing equations vere modified for the . program to 

permit direct solution for A by substituting: 

i 

A • A ^ 

B . z . ds 

1 1 


The sets of gas and liquid phase equations are coupled through mass and nK>men- 
tum exchange between phases. For droplet gasification, the simple evaporation 
coefficient model is utilized: 



REPRODUCIBITjr: OF THIS 
ORIGINAL PAGE IS POOR 


where the evaporation coefficient is 



and 

Nu.” = 2 + 0.53 Re " 
J J 


Drag forces on spray droplets are expressed by 





Pi 



( 23 ) 


( 24 ) 


( 25 ) 
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with the drag coefficient specified as 

C- " = 2'i (Uc .") : lie s 80 

J J 

0.217 

- 0.271 (no.") : Re." > 80 

J J 


( 26 ) 


2 . 2 . 2^1 Performance Parameters 

Two separate parameters are calculated which are indicative of the overall 

degree of propellant mixing* These are calculated once in LISP, based on the 

flowrates associated with the LISP mesh points, and once in STC, based on the 

initl'^l flowrates supplied to the stream tubes. One parameter is E , a mixing 

m 

efficiency factor due to Rupe (Ref. 7) which expresses a mass-weighted average 
approach of local oxidizer mass fractions to the overall injected mass fraction: 


E 

m 


100 


1- E 


Wj^(R-r^) 

Wr 


n 


E 

1-1 


Wi(R-ri) 

tJ(R-l) 


where : 

n * number of samples with R^r 
n * number of samples with R<r 
w = local propellant flowrate, Ib^/sec 
ft * total propellant flowrate, Ib^/sec 
r,r = local oxidizer mass fractions, w^/w 
R “ injection oxidizer mass fraction, W^/ft 


(27) 


The second parameter is a mixing c* efficiency, which represits the 

maximum attainable c* efficiency corresponding to complete propellant gasification: 


n+n 

c*(c )w, 

* 1=1 

c*(c 

Inj 


(28) 


where is the injection mixture ratio (W^/Wp) , is local mixture ratio 

(wo^/wpj^) and w^^, ft, n and n have the same meanings as above. TheoreMcal 
characteristic velocity is tabulated as a function of rixture ratio. 
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During STC’s multiple stream tube analysis, a single value of efficiency 
is calculated from the n^ stream tubes’ data at the throat plane: 


n 


c* 


c*(Cinj)W 


( 29 ) 


Note the distinction between Eq. 29 and 28: local gasified propellant mixture 
ratios and flowrates are used in Eq. 29 rather than local total mixture ratios 
and flowrates, as in Eq. 28. 


2.2.3 Method of Solution 

The numerical integration scheme used to solve each stream tube’s system of 
equations is the simplest first-order Runge-Kutta (or Euler) method. Selected 
for its simplicity, minimal data storage requirements, low execution times and 
numerical stability, this method’s accuracy is strongly dependent upon using 
sufficiently small step sizes. This limitation is reduced in importance by 
using backwards differencing ^n writing finite-difference equations and by 
solving the equations twice, using predicted values from the first, or predictor, 
solution as input data for a second, corrector, solution. 

The STC program is first run in a single stream tube mode, i.e., a one-dimen- 
sional subsonic combustion analysis is made for the entire chamber using 
appropriate sums and averages of initial stream tube variables. This is done 
for two reasons: (1) to verify consistency of input data (initial-plane pres- 
sure is adjusted until the cne-dimensional throat velocity is within a small 
tolerance of the calculated throat sound speed) , and (2) to provide a mean 
adiabatic expansion coefficient, y# for combustion gas flow in the convergent 
part of the exhaust nozzle. 


The latter coef ficien t_is gi ven by: 


\ pi- 


(30) 


where the subscript 1 refers to the beginning of nozzle convergence, the 
variables p* and p* are at sonic conditions and the over-bars refer to the 



one-dimensional flow analysis* It is used by the TRANS computer program 
(described in the next subsection) to calculate the coordinates of constant 
pressure surfaces (isobars) for transonic flow in the nozzle* TRANS isobars 
are generated and transferred to STC in non-dimensional terms, so their use in 
STC requires knowledge of the nozzle throat radius, Rj, (an input parameter) , 
and sonic flow pressure, p*. Discussion of the evaluation of p* is included 
in the following outline of STC*s multiple stream tube solution. 

Following STC average single stream tube analysis and TRANS analysis, the 
initial plane is reinitialized with its original input and the STC program 
is run in a multiple stream tube mode. Sequentially: 


1. The main iteration loop performs the z direction marching. It begins with 
estimates that the changes in chamber pressure, stream tube gas velocities 
and densities across the next Az increment are equal to their gradients 

in the preceding single stream tube analysis* 

2. These estimated properties are used to calculate predicted values for all 
spray behavior parameters. Drag, evaporation and other spray droplet 
size group parameters are computed with controls to: (a) limit evaporation 
if it is found to exceed the amount of spray available, and (b) avoid 


having the sign of 


u ~ (Uj ”) 
s s 


change in a non-physical way due to 


. 


overestimation of drag forces, i*e., droplets cannot accelerate or 
decelerate past the gas velocity within a given Az* 

Evaporated spray weights are added to the previous gasified propellant sums 
and the gas phase mixture ratios are computed* Combustion gas properties 
corresponding to these mixture ratios are obtained by linear interpolation 
in the properties table. Corrected estimates are then made of gas 
temperatures^ densities, and stream tube areas. 

Spray gasification and draf, terms are next treated as known constants in 
an implicit-explicit two-step solution for the gas phase properties, stream 
tube flow areas and chamber pressure in plane Z 2 « In the first step, a 
pressure level in plane z^ is assumed in order to calculate a predicted 
distribution of stream tube areas in that plane* In the second step, that 
distribution of areas is assumed to be valid, making possible an explicit 
solution for the gas temperatures, densities and pressure at z^. 
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5. The foregoing paragraphs 2, 3, and 4 describe a predictor cycle. Their 
calculations aie repeated in a corrector cycle, using the predictor cycle’s 
calculated z^-plane results instead of the estimated properties (paragraph 1) . 
If evaporation of a spray group is calculated in the corrector cycle to 
exceed the total mass of that group, the group is said to be completely 
gasified. (For improved accuracy, a user may elect to perform additional 
corrector cycles.) 

6. At this point, normal program flow Is printout of computed data at plane 

reinitialization of plane z^ at pl^nc z^ for a new Az and progressing 
again through the paragraph 1 through paragraph 5 computations. This 
procedure would continue for non-axisymmetric stream tubes until nozzle 
throat plane was reached. With DER’s axisymmetric annular stream tubes, 
however, the procedure is changed (as follows) at tiie position in the 
nozzle where curvature of the isobaric surfaces is introduced, and marching 
continues for as many as 25 Az’s past the throat plane. 

7. An approximate value of p* is estimated from the nozzle threat plane pres* 
sure of the preceding averaged, single stream tube analysis: 

p* a p* p(z^/p(z^). By multiplying the reduced pressures, p/p*, of the 
TRANS isobars by this value of p*, absolute pressures are calculated for 
the transonic flow field. These are imposed upon the multiple stream 
tube nozzle flow. 

The furthest upstream TRANS isobar may be planar or curved, depending upon 
the nozzle’s radius ratio and shape of its convergent section. If it is 
curved, it is desirable to introduce a gradual transition from planar 
isobars to that first curved isobar which the solution encounters. Also, 
a gradual transition is desirable to smooth out any discontinuity in 
pressure levels between those solved for upstream and those imposed down- 
stream of the transition. The gradual transition is provided by stopping 
the solution for pressure level at a position that is upstream of the 
nozzle throat by 1.3 times the axial distance that the furthest-upstream 
TRANS isobar intersects the nozzle wall, aad using linear interpolation 
to impose absolute pressures over the transition interval. 
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In the transonic region, imposition of absolute pressures makes the 
paragraph 4 solution tor pressure level redundant, sc that the gas phase 
calculation is reduced to a one-step, explicit solution for stream t be 
areas. Because absolute pressures have been imposed, the solution rov 
provides absolute values of stream tube areas and those may or may not . m 
to the local nozzle flow area (i.e., satisfy area continuity). By this 
method, area continuity can be satisfied only by finding the appropri ite 
value of p* to define the proper nozzle pressure level. This is done nly 
for the minimum flow area, as follows. 

8. As the solution marches through the transonic flow regime, the minimum 
value of the sum of the calculated stream tube areas at any preceding 
z-plane is stored. Later, after the solution has reached a z-plane that’s 
wholly downstream of the TRANS isobar which intersects the nozzle wall at 
the throat, this minimum area sum is compared with the input geometric 
area of the nozzle throat. It is the match of these areas that constitutes 


satisfaction of the nozzle throat boundary condition. If the fractional 
deviation is less than some input tolerance, c , the solution is complete 
and data for subsequent input to the TDK computer program are punched out. 


If the deviation lies between one and three timesGA , the value of p* is 

t 

adjusted and the solution is recalculated from the point in the nozzle 


where absolute pressures were imposed. For deviatiens exceeding three ^A^* 
the STC Initial-plane pressure is adjusted and the entire multiple stream 


tube solution is recalculated. 


2.3 TRANSONIC NOZZLE FLOW: 

TRANS SUBPROGRAM BLOCK 

A transonic flow analysis section was adapted from the reference TDK computer 
program (Ref. 8), as modified (Ref. 9) to utilize an elliptic coordinate trans- 
formation solution method (Ref. 10). This section was removed from the TDK 
program and modified so that it would generate a family of isobaric lines 
throughout the transonic flow regime and provide a computer-plotted graph of 
that family. The necessary input data are obtained from the averaged, single 
stream tube solution of STC, so this TRANS subprogram block gives a homogeneous 
flow solution* For homogeneous flow, TRANS solutions are stable with radius 
ratios as small as 5/8. 
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and a 


As input data, TRANS needs values only of the nozzle throat radius, 
mean expansion coefficient, y* Isobaric coordinates are calcvlated in terras of 
axial distance, X, from the throat plane and radial distance, R, from the 
nozzle axis; both dimensions are normalized to the throat radius, haltin^e 
isobars are generated, one at a time, by starting downstream of the throat and 
marching upstream with equal intervals, Aa, in the angle between the nozzle 
axis and a line tangent to the nozzle wall at the isobar/wall intersection 
point. The program is structured such that that intersection point for the 
fifth isobai Is at the throat; this isobar later becomes the TDK sta^.^li. . 
Four isobar/wall intersection pc .nts He downstream of the throat (a>0) and 
the remainder lie upstream of the throat (a<0). lac anf;ular im, ,oi between 
isobars is given by: 

Aa . - 

Generation of isobars continues until either (1) there are twenty of them or 
(2) an isobar exhibits significant reverse, or upstream curvature. In the 
latter case, that last upstream-^curving isobar is replaced with a planar 
surface. 

Two computer-plotted examples from TRANS analyses are shown in Fig. 3, where: 
the nozzle axis is at the bottom (R/RT =0); flow direction is from left to 
right; a portion of the nozzle wall, defined by -JO® - wall angle 5 + 8®, is 
shown as the upper curve; isoLars are generated from right to left at nozzle 
wall angle intervals of 2.00® (Eq. 31 was not used in these runs). The 
monotonic downstream curvature of the constant nressure surfaces is apparent, 
as is its accentuation by lowering the nozzle radius ratio. Included on the 
figure are tables which list the pressure ratio, p/p*, and Mach number for 
each isobar. 

The TRAILS program also calculates a nozzle discharge coefficient using the 
3rd order equation given in Ref. 10: 


(31) 
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Before proceeding to multiple stream tube STC program analysis, the STC initial- 
plane pressure which satisfied one-dimensional sonic throat flow is divided 
by to obtain an improved estimate of initial plane pressure and expedite 
convergence on the nozzle throat boundary condition. 

2.4 SUPERSONIC NOZZLE FLOW: 

INTERFACE WITH THE jm COMPUTER PROGRAM 

The previous version (Ref. 2) of DER included a TIK subprogram block for per- 
forming supersonic nozzle expansion of multiple, axisymmetric , gaseous stream 
tube flows. That subprogram block was based on minimum modifications of the 
TDK long-form option, described in Ref. 8, wherein only supersonic analysis 
is performed. A supersonic start line was initialized from STC computed data 
and from equilibrium computations using the ODE (one-dimensional equilibrium) 
section of TDK. The TDK start line data generated by STC were printed out at 
the end of STC computations and a corresponding punched-card deck was generated. 
Thus, TDK could be run, optionally, either in sequence with STC using data 
transferred via scratch data units or separately using punched card input data. 

A substantial discrepancy was found to exist between the start line gas 
temperatures and densities computed by STC and those computed by TDK as start 
line equilibrium. The discrepancy is caused by the difference between STC’s 
frozen expansion from stagnation equilibrium and TDK’s static equilibrium. 

Under another program (Ref. 11) a modified STC solution method has been 
developed which eliminated this discrepancy. This method involved: expanding 
the input tables of equilibi^ium gas properties to be functions of both mixture 
ratio and flow Mach number; using static equilibrium properties rather than 
stagnation; forming (internally) pseudo-stagnation temperature arrays from 
the input static temperatures; and, interpolating to local mixture ratio and 
Mach number in these tables. 

At the same time, however, an improved version of TDK was also developed (Ref. 3) 
which is capable of performing kinetic expansion analyses for multiple gaseous 
stream tubes In subsonic and transonic flow regimes as well as in the super- 
sonic nozzle flow. While it would have been possible to adopt the modified 
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STC subprogram block from Ref. 11 and modify the improved TDK to accept STC 
data as input, these actions were not taken. Rather, the TDK subprogram block 
was removed from DER and the frozen expansion from stagnation equilibrium was 
retained. A subroutine was added which prints and punclies, in NAMELIST format, 
the throat plane data needed from STC to continue the multiple stream tube 
analysis via the improved TDK computer program. 


Throat plane data punched out (and the TDK parameters to which they correspond) 
are: the number of stream tubes (NZ0NES) , a stream-tube-area-weighted mean 
stagnation pressure (P(l)), 


P * 

o 


1?1 'o** 

nq^ 

i«l 


(33) 


and, for each stream tube, the gasified propellant mixture ratio (0FSKED) and 
mass fraction of the total gasified propellant flow within that stream tube (XM). 


It is anticipated tha^. these stream tube data may be used to initiate TI^ 
nozzle expansion analysis at the throat plane or, if kinetic effects are 
believed to be important in the subsonic flow, at some plane upstream of the 
nozzle throat. Thus, it is now possible to overlap the spray gasification 
analysis of STC with the non-equilibrium combustion analysis of TDK, but in an 
uncoupled manner. The location of the TDK initial plane is designated by the 
program user, by specifying the contraction ratio (ECRAT) for that start plane 
and also the subsonic area ratio for that plane (SUBAR(l), equal to the 
contraction ratio). This non-physical stream tube data transfer is illustrated 
in Fig. 4. Other data needed for initializing TDK are*: reactant cards (pages 
6-9 and 10), reactions cards (pages 6-27 to 31), nozzle design parameters 
(pages 6-32 to 34), integration and print controls (pages 6-35 to 37) and 
TDK controls and nozzle divergence geometry ipages 6-42 to 46) . 


*Page numbers in parentheses refer to the Program User’s Manual Section of 

Ref. 3 . 
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Modifications have been made to the improved TDK program to make it compatible 
with DER. TDK analyzes only the gasified propellant flow. Account is taken 
of the propellant mass loss represented by residual sprays passing through the 
throat, assuming no continued evaporation or acceleration that such sprays 
might undergo in the supersonic nozzle expansion section. Residual spray is 
not added to the TDK initial plane momentum. 
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3, DER PROGRAM USER'S GUIDE 


The DER computer program was developed for operation on Rocketdwe’s IBM System 
360, Mod 50/65 computer which is designed to run programs written in Fortran H. 
So that this program would be compatible with other computers, however, it was 
written in Fortran IV (which is a subset of Fortran H) . There are, of course, 
some sub-programs which may not be operable on other than the Rocketdyne com- 
puter; for most other computers, these are probably restricted to the data- 
plotting functions and can be replaced by dummy subroutines without detriment 
to the rest of the program functions. In order to run the program on any 
computer, a user must supply program control cards that are compatible with 
his compiler, link editor, etc. The program makes extensive use of overlay in 
order to reduce computer storage requirements; the overlay structure is 
described in Section 3.1. Storage required and other operational considerations 
are discussed in Section 3.4. 

Operation of the DER computer program also depends upon a user-supplied data 
deck, through which he specifies details for the particular combustor and 
propellants he desires to analyze. Assembling a data deck usually involves 
expenditure of a substantial effort. There are separate sections of the data 
deck for each major subprogram block; they are assembled in the order in which 
the program calculations proceed. Details concerning data deck assemblage are 
given in Section 3.3. 

3.1 PROGRAM OVERLAY STRUCTURE 

The DER computer program can be overlayed to operate on a computer having only 
moderate storage capacity. The recommended overlay structure is shown in 
Fig. 5, where the branches are set us such that most of them are required only 
once. The program length with this structure is 42,450 words; maximum length 
occurs when the branch comprising CSPRAY, et seq., is loaded. Total storage 
used during execution on Rocketdyne*s IBM System 360, Mod. 65 computer is 
about 47,200 words. 
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ORIGIN A 


LISP 

/ELEM/ 

/SHAPEC/ 

/LABML/ 

/CPLJT/ 
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ORIGIN B 
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PL0TN 
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Figure 5 


ORIGIN C 

CINPUT 

STAPE 

SCRMBL 

S0RTR4 

/SCR/ 

/I SCR/ 

WVD 


ORIGIN D 


im 

TRANS 

KPRIME 

N4NAIN 

ESUMB 

HALL 


ITER 


PL0H 

1 

CDTRAN 


CSPRAY 

EVAPS 

CGDYN 

QUAD 

CPRINT 

AV030A 

^TER8 

PVSST 

PVSR 

PUNAME 

/PUNNAM/ 


DER Prograza Overlay Structure 
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Label C<8 mM 0N blocks needed only within a particular segment are included in 
the overlay structure# Each is denoted in Fig# 5 by enclosing the C0MM0N 
block nanie between /'s, e.g. ,/CNTRE/. There are several label C0MM0N blocks 
which are not shown in Fig. 5; these are understood to be part of the root 
segment, l.e., that including DER, HEADER, et seq. 

3.2 PROGRAM SUBROUTINES 

Brief statements are given about each of the subroutines in the DER computer 
program, indicating the major functions performed# These are intended to assist 
the interested user in finding quickly where particular operations are per- 
formed. No attempt is made to give detailed descriptions or to delineate 
fully the equations solved or the solution methods used# The subprograms are 
divided into three groups, corresponding to the main stem (DEF , et seq.), the 
LISP branch and the STC branch, respectively, of the overlay structure. The 
sequence of subroutines in this section conforms to the structure of the 
overlay chart. Fig. 5. 

3.2.1 DER (Main Stem) Subroutines 

The DER main program is a brief executive program that calls subroutine HEADER 
and, based on values read-in there for certain control parameters, calls 
sequentially subprograms LISP, STC, TRANS and PUNAME. If a nonzero value of 
the integer INERR is returned from LISP, the rest of the calls are bypassed 
and execution is terminated. 

Subroutine HEADER reads a set of four comment cards, descriptive of the 
particular case being run, and a set of program flow control integers. It 
then prints out a title page for the case being run. 

Subroutine TIME prints out the date, time of day and elapsed clock (real) time 
since the previous call of TIME. This subroutine calls library subroutines 
CL0CK and CDATEV; if either of these Is unavailable in the library of a 
particular computer, subroutine TIME should be replaced by one which utilizes 
available library routines or by a dummy subroutine TIME. 
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Subroutine SCALE determines rv**quired limits and increments for scaling the 
axes of computer^plotted graphs. It is referred to by subprograms CLEVLS, 
LGRID, PL^TT, and LLGRID. 

Subroutine LINEG connects successive points of a computer-plotted function 
with straight line segments. Library subroutine L3,;\i:.V is called by IINEG. 

Function subroutine Y0F (X, XT, YT, N, NP) performs an NP-point Lagrangian 
interpolation to determine Y«Y0F ■ f(X) corresponding to the point (X,Y) in 
the N-point array (XT, YT). Subroutine L0CATE is called by Y0F to find the 
first of the points (XT(K), YT(K)) for interpolation, where 1^ K%. Sub- 
routine XITRP is called by Y0F to perform the actual point-to-point inter- 
polation. 

Subroutine ITRP2 (XI, Tl, Nl, II, X2, T2, N2, 12, YT, LI, Y) performs a 
double Interpolation to find Y«f(Xl,X2) from the YT(N1,N2) array corresponding 
to the position XI in the Tl(l) to Tl(hl) array and the position X2 in the 
T2(l) to T2(N2) array. Subroutine L0CATE is called to determine the appro- 
priate first points and subroutine XITRP performs the actual interpolation. 

3.2.2 LISP Subroutines 

Subroutine LISP is quite large and, as the executive subroutine of this sub- 
program block, performs many functions. The first of these is reading-in 
all of the data required for any of its associated subroutines. Some of the 
input data are tested, as they are read-in, for consistency. If there are 
inconsistencies in the data which would invalidate LISP*s computations, an 
error message is written out and the case is deleted immediately. Given 
acceptable data, LISP next sets up the coordinates of the network of mesh 
points in plane Zq. Then, sequentially, subroutines EFL0W, DSIZE, SC0EFF, 
CCKP, FANIN, MFLUX and BNDY are called to perform their functions, if appro- 
priate. Subroutine LISP then scales the computed mesh point mass fluxes to 
ensure conservation of total propellant flowrates, writes out the computed 
data (both before and after evaporation), generates a scratch data record of 
information to be transferred to STC, calls subroutine MIXSFF, calculates 
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and prints mean properties for a combustion gas flow at Zq resulting from gasi- 
fication of propellants upstream of that plane. Subroutines PL0TN, CRTLBL, 
PL0T2 and PL0TC are called ' f controls have been set to generate computer 
plotted graphs of data and, finally, the entire LISP procedure is repeated if 
analysis at a second (or third) LISP collection plane has boen called for by 
reading in non-zero Z0M2 (and Z0M3) . 

Subroutine EFL0W calculates propellant injection flowrates and velocities for 
each injection element of Types 1 through 5 and 7 and 8. For the gas/liquid 
element types (Types 6, 9 and 10), the injection flowrates are supplied as 
input data. Variables denoting the flowrates of each propellant from each 
orifice (or equivalent orifice) of each element are also assigned. 

Another function performed in EFL0W is calculation of correlation coefficient 
a for the gaseous (fuel) propellant of a Type 6 (gas/liquid coaxial jet) 
element. This is done on the basis that the vaporized portion of the oxidizer 
spray is diffused and mixed uniformly with the gaseous fuel and that gas phase 
momentum is conserved. 

Subroutine DSIZE contains equations for calculation of propellant mass median 
droplet diameters for elements of Type 1 through 5. These compotations are 
bypassed if non-zero droplet diameters are supplied as input drta. 

Subroutine SC0EF provides, for Type 1 through 6 ana Type 9 elements, val^ 's 

of the mass flux correlation coeffi< Lents, C. through C . , a and b for both 

1 o 

propellants. These variables are read-in for Type 8 and 10 elements. SC0EF 
calls upon appropriate coefficient library subroutines FSA through FSAC0L. 

Subroutine CCKT calculates modified evaporation coefficients, (CKPl and 2), 
in the event that these variables were both read-in with values greater than 1. 
In that case, the read-ln d ita are taken to be the propellant latent heats of 
vaporization and this subroutine replaces them with a calculated . If 
both read-in values are less than unity, CCKP is bypassed. 
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Function subroutine FSA provides values for correlation coefficients a and C. 

4 

for Type 1 (unlike doublet) and 2 (like doublet) elements. Function sub- 
routine FSB provides values for coefficients b, C, , C_ and C, for Type 

L A j j o 

1 and 2 elements. Both propellants are included. 

Function subroutines FPDSA and FPSSA provide values for coefficients a and 
for the outer orifices and the central orifice, respectively, of Type 5 
(liquid/liquid 4-on-l) elements. 

Similarly, function subroutines FTDSA and FTSSA evaluate coefficients a and 

for the outer (fuel) orifices and the central (oxidizer) orifice, res: ,tively, 

of Type 3 (liquid/liquid triplet) elements. In like manner, coefficients b, 

and C, for the triplet *s outer and central orifices are obtained from 
A o 

function subroutines FTDSB and FTssB, respectively. 

Spray distribution coefficients for Type 3 (like doublet pair) injection 
elements are obtained from function subroutines FLDPSA, FLDPSB and FLDP0A. 

The a coefficients for each propellant differ in the (-x) direction from those 
in the (+x) direction. The value of a for the "same side” of the origin as a 
given propellant’s doublet lies is obtained from FLDPSA and that for the 
•opposite side” of the origin from FLDP0A. The b coefficients are calculated 
by FLDPSB. 

Function subroutine FSAC0L calculates coefficients a and for the central 
liouid jet of a gas/liquid coaxial jet element (Type 6). 

Subroutine FANIN was originally intended to account for interference between 
the spray trajector'^es of closely-spaced, neighboring spray fans. This 
function has not been "»rcvided, however, so the subroutine’s only function 
is assigning and writing out values of coefficient a foi positions to the 
right (x>0) and left (x<0) of the element origins. 

Subrov^tine MFLUX performs primary functions of calculating each injection 
element’s contributions to each collection plane mesh point, summing them 
to obtain mesh point propellant fluxes (both before and after accounting for 
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partial spray gasification) and calculating mean mesh point spray velocity 
vectors. Secondary functions are concerned with computing and printing some 
variables that are indicative of spray flowrate collection efficiencies. 

MFLUX calls upon a number of computer library trigonometric functions. 
Arithmetic statement functions define functions SIND(X) and C0SD(X) for 
angular arguments X in degrees; these definitions should be removed if the 
library contains these functions. 

Subroutine BNDY provides for folding propellant flows \diich fail outside of 
the defined chamber segment being analyzed back into particular mesh points. 
Flow computed to pass radially outward through the chamber wall is folded 
back into the wall mesh points. Treatment of flow which passes azimuthally 
through radial planes of symmetry depends upon the value of JSYM and upon 
whether there are physical barriers (e.g., radial baffles) there. BNDY calls 
upon subroutines SUMM to sum the propellants folded into the outer wall mesh 
points, SUH2 to accomplish the folding at radial lines of symmetry and SUMV 
to calculate appropriate adjustments of the propellant mean velocity vector^ 
at mesh points which have received folded flows. 

Subroutine MIXEFF computes values of Rupees mixing efficiency factor, E , and 

" m 

a mass weighted mean characteristic velocity efficiency, mix* based on 
mesh point propellant flows and a read~in table of c’*' vs mixture ratio (0/F). 

For Type 6 or 10 gas/liquid injection elements, subroutine GFL0W calculates 
and prints out local mesh point combustion gas conditions. It is assumed that 
the vaporized portions of the liquid sprays are mixed with the initial gaseous 
propellant at uniform mixture ratio but the gaseous mass flux can be non- 
uniform, Combustion gas properties are obtained as functions of that mixture 
ratio by interpolation in read-in tables. 

A package of subrout !nes is required for punching data cards for 1f)C. LISP 
calls PUN to produce cards to specify spray drop size, velocity components and 
flowrates at each mesh point. Also, it calls GPUN to produce cards to specify 
gas mixture ratio, density, velocity and temperature at each mesh point. 
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Subroutine PUN uses some qeneralized punch subroutines: PPUN» XPUN, XPUNI and 

XPIJN?. A variable format is used in PPUN which is set up by subroutine S4MAT. 
This latter subroutine Is dependent on the type of computer used, i.e., in 
respect to the number of Hollerith character^ stored in each word location. 

'I’hc remain! HR LISP subroutines arc all concerned with generation of computer- 
plotted graphical output. They are designed to utilise Stromberg-Carlson 
SC-4010 equipment in producing 9x9-in CRT plots. Subroutine P1.0TN produces a 
cross-scct ional graph of the chanfcer section being analyzed if NCRT ^ 0. The 
location of each element ‘s origin is indicated thereon. Subroutine PL0T2 


(Continued on next page) 
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1b called if NCRT > 0;it plots propellant mass fluu vs 0 along specified 
r = constant arcs* The basic rectangular grid Is established by subroutine 
LGRIf). Subroutine CRTl.BL writes labels on the abscissa and ordinate of each 
PL0TN or P1.0T2 graph* 

Subroutine PL0TC constructs circular plots of the system of LISP mesh points 
and» on success I Vi* graphs » plots contour level maps of fuel» oxidizer and 
total propellant flux and one of ci function of the basic mixture ratio* If the 
maximum and minimum magnitudes of the contour variables are not specified, 
they are determined by subroutine MXMN* Subroutine CLEVLS establishes the 
actual contour levels plotted, assigning some convenient interval between 
levels* Subroutine C0NTUR is called to do the actual contour line plotting. 

It, in turn, calls upon subroutine C0NTRE to establish the contours position 
between four mesh points and C0NTRP to plot and label the contours in 
repeating pic-slice-shaped chamber segments* 

3*2*3 STC Subroutines 

Subroutine STC is the main executive subroutine for the STC subprogram block 
of DER* After calling subroutines CINPUT and CINIT for reading-in data and 
initializing the start-plane conditions, respectively, a main D0-loop is 
established which performs the marching from z^ through the nozzle throat* 

Within this D0, the combustion model equations are solved by stepping from 
\nown conditions in one plane and calculating conditions in a second plane 
a small Az downstream. A number of other subroutines are called upon to carry 
out the actual computations* Flow control for predictor-corrector cycle 
calculations resides in STC. The logic for testing to see if the single 
stream tube analysis has converged on a solution, for reinitializing multiple 
stream tube conditions and for testing whether the multiple stream tube analysis 
has converged are also in STC* 

Subroutine CINPUT reads all of the punched card data needed as input by the 
STC subprogram block and, if appropriate, calls upon subroutine STAPE to 
read the scratch data unit recorded by subroutine LISP and set up the initial 
stream tube areas, propellant flowrates, gas phase mixture ratio and spray 
properties. If the option Is specified for starting the STC analysis from 
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3DC data, STAPE reads data from cards punched by 3DC rather than from a scratch 
data unit and droplet flowrate, velocity and diameter are averaged at each mesh 
point by subroutine WVD. 
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Subroutine STATE reads data transferred from LISP via scratch data unit 2 and 
uses it to set up the stream tubes* initial cross-sectional areas and propel- 
lant flowrates and conditions at the SIC start plane. STATE calls subroutine 
SCRHBL to perform the assignment of mesh points to particular stream tubes 
and to do the appropriate summation and averaging of areas and propellant 
parameters. Data returned from SCRMBL are multiplied by an integer scale 
factor to convert from a pie-slice-shaped sector to a full circular chamber 
cross-section. Stream tube initialization data are punched-out so they may 
be used as input for subsequent STC subprogram block runs. 

Subroutine SCRMBL is a fairly complicated subroutine, as it contains the logic 
for accomplishing the several optional ways of combining mesh point flows into 
stream tube flows described on page 8 . When mesh points are to be combined 
on the basis of mixture ratio, subroutine S0RTR4 is called to sort the mesh 
points into a sequence with increasing mixture ratio. This is done by tagging 
each mesh point with one value of a subscripted integer label IC(K) whose 
subscript K is uniquely related to the mesh points r,9 subscripts. 

Subroutine CINIT performs several functions in preparing for STC’s main 
D0-loop marching solution. Dummy arrays are defined for saving stream tube 
initialization data for subsequent occasions when the start plane is reini- 
tialized and STC*s computations are restarted. Subroutine KTRIME is called to 
compute, store and print tables of propellant evaporation coefficients as 
functions of combustion gas static temperature. Subroutine AVAR is called 
to set up a table of chamber cross-sectional area vs chamber length and to 
define some other geometric parameters* If STC single stream tube analysis 
is being prepared for, the stream tubes are (temporarily) combined into one; 
thus, only the first of NST stream tubes is properly defined for this option 
and it carries the total injected propellant flowrates and occupies the entire 
combustion chamber cross-section. CINIT also contains the logic for re- 
establishing the multiple stream tube conditions upon completion of single 
stream tube analysis. Other functions include initialization of the stream 
tube combustion gas properties and printing of some initial stream tube data. 
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Subroutine TABLES interpolates in read-in tables of combustion gas properties 
at local stream tube mixture ratio to find local combustion gas stagnation 
temperature, gamma, molecular weight and viscosity. Local stagnation sound 
speed and throat static sound speeds are calculated. 

Subroutine RC00RD establishes the radial coordinates of the dividing stream 
lines betx^een stream tubes and the incremental stream tube path length cor- 
responding to a Az axial increment. Downstream of the nozzle throat, the 
dividing stream lines* intersections with the fifth TRANS isobar are solved for 
and printed out. 

Subroutine STCRT records the dividing stream line radii as functions of chamber 
length and, upon completion of the chamber length marching analysis, generates 
an axial chamber cross-section plot. 

Subroutine CHANGS reassigns the values of stream tube parameters in the 
upstream z-plane as being equal to those in the downstream z-plane, as prepara- 
tion for another Az step* (When called from CINIT, subroutine CHANGS* function 
is inverted, assigning upstream conditions to the downstream plane.) 

Subroutines KPRIME and AVAR were commented upon along with subroutine CINIT. 

Subroutine ESUBM calculates a value of the Rupe mixing efficiency factor, E , 

^ ' m 

for the stream-tube-striated combustor flow. Total stream tube propellant 
flows are considered, so this value corresponds to complete propellant 
gasification. 

Subroutine TRANS is the executive program for the TRANS subprogram block. It 
is called by the DER main program. TRANS defines a number of control 
parameters, then writes out a table of input and control data before calling 
subroutine N4MAIN, which perfoims the transonic analysis. 

Subroutine N4MAIN solves for the radial and axial coordinates of points on 
isobarlc surfaces in an axisymmetric nozzle’s transonic flow regime. Sub- 
routine Hall is called upon to perform a modified Hall solution of transonic 
flow equations for a homogeneous constant flowrate gas. Subroutine ITER is 
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called to solve for the axial coordinate roots. The isobar io surfaces’ 
coordinates are plotted by subroutine PL0TT and are tabulated, together with 
each line’s reduced pressure, Mach number and flow directions, by N4MAIN. 

Subroutine CDTRAN is called by TRANS to calculate and print a nozzle discharge 
coefficient for the foregoing transonic solution. 

Subroutine CSPRAY calculates the behavior of propellant spray si::e groups. 
Subroutine EVAPS is called to compute the size group gasification rates. 

CSPRAY then determines whether those rates are permissible, adiusts them if 
necessary, reduces the liquid spray flowrates and calculates the size group 
droplet diameters and velocities. If there are multiple size groups of a 
propellant, the droplet number flowrate is constant and droplet diameter is 
a variable. However, if there is a single size group for a propellant, the 
diameter is held constant and the number flowrate is reduced as gasification 
proceeds* 

Subroutine EVAPS calculates size group gasification. In doing so, it calcu*- 
lates droplet Reynolds numbers based on mean droplet film properties. Film 
viscosity is approximated via Wilke’s equation. 

Subroutine CGDYN solves the stream tube gas dynamic equations. Four finite- 
difference equations (continuity, momentum, state and adiabatic) in four 
variables (velocity, temperature, density and either area or pressure) are 
combined to form a single quadratic equation for gas velocity* A gas velocity 
root is determined explicitly, via function QUAD, and is used in the individual 
equations to solve explicitly for the other variables. Different quadratic 
equations are used, depending upon whether pressure or stream tube area is 
(assumed to be) known at plane 2, In the multiple stream rube case, pressure 
is first assumed to be known in plane 2 in order to solve for a distribution 
of areas. In a second step, that distribution of downstream areas is assumed 
to be valid and the downstream pressures are solved for; an area-weighted 
mean pressure is assigned to the plane. Obviously, in the case of a single 
stream tube, the downstream area is known and the first step is not needed. 
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similarly, imposition of absolute pressures, based on the TRANS distribution 
in the nozzle, makes it possible to bypass the second step for the multiple 
stream tube case. 

2 

Function subprogram QUAD (A,B,C,X1) solves the quadratic equation AX + BX + C*=0 

for the root whose value is closest to the argument Xl* If the discriminant 
2 

B^-4AC < 0, it is set equal to zero and the single root is returned. 

Subroutine CPRINT’s function Is to print computed stream tube pressure, area, 
gas stream and propellant spray data at specified multiples of Az from the 
STC start plane. Different formats are used for the single and multiple 
stream tube cases. 

Subroutine AVD30A calculates a value of volume-number mean droplet diameter, 

Djq, for each propellant. It is called by subroutine CPRINT for single stream 
tube print out. 

Subroutine ITERS records data, on scratch data unit No. 3, at the plane 
immediately preceding that plane where a radial distribution of pressures is 
imposed on the multiple strea*** tube flow. These data may be used to re- 
initialize stream tube combustion analysis with adjusted values of pressure, 
gas densities and velocities at that plane if the computed minimum stream tube 
area sum is not sufficiently close to the nozzle throat area. Another data 
record is recorded at the throat plane; it is read by subroutine PUNAME and 
used for preparing input data for the improved TDK. 

Subroutine PVSST calculates the radial distribution of pressure in a z-plane 
and assigns a pressure to each stream tube corresponding to its mean radial 
position. Subroutine PVSR is called to perform the interpolation between 
TRANS isobars. 

Subroutine PUNAME reads data from data unit No. 3 and calculates data for 
initializing the Improved version of TDK. A punched card deck is generated, 
using the NAMELIST format appropriate for TDK input data. The deck is also 
printed out. Additionally, values of a number of other variables, which may 
by needed in a future revision of TDK, are calculated and printed out. 
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3.3 PROGRAM INPUT DATA 


In general » each of the two major program blocks LISP and STC performs its own 
input and output functions, and they can be run sequentially or individually. 
When executed sequentially, coupling data are transmitted via scratch data 
unit from LISP to STC and other data are supplied directly to STC as required 
via punched card input. When run individually, all data to a block are supplied 
via card input. Punched card output of the coupling data are also provided 
(optionally for STC) for convenience in running subsequent STC analyses 
without having to rerun LISP and for subsequent analyses using the improved 
TDK program (Ref. 3). (For example, the same LISP output might be used for 
several STC runs v^ith different chamber and nozzle lengths, or TDK might be 
run more than once with different nozzle contours or expansion ratios using 
a single STC output.) 

Input data requirements are indicated in Tables 1 and 2. Table 1 shows, in 
the required order, Fortran code sheets for all of the required data cards. 

The ’’DESCRIPTION’* column gives the DER program’s Fortran name for each 
variable; definitions of these input variables and their units are listed 
alphabetically in Table 2. Values of the variables are entered in appropriate 
fields occupying columns 1 through 72, using the format designated in the 
’’IDENTIFICATION” space at the bottom of each card in Table 1. Sequence numbers 
are denoted in columns 73 through 80 which both leave room for the maximum 
sizes of subscripted arrays and ensure that the data deck is sortable. In 
some instances, where multiple sets of the same cards may need to be entered, 
blank spaces appear in the sequence numbers. 

Whether or not certain data are actually required for a particular case depends 
upon values read-in earlier in the case particular control parameters. 

These controls are Indicated in Table 1 either as a condition attending major 
blocks of cards or parenthetically beside the card specification. 

Discussions guiding the selection of input data are subdivided under LISP 
and STC subheadings following Table 2. 
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TABLE 1. DER Computer Program, Input Data Cards 


Data Read by Subroutine Header, Every Case 



3 



















TABLE 1 (COrJT'D) 




(if ILISP = 0 
omit cards 



2010 through 
4630) 













TABLE 1 (CONT'D) 
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AO 



























































































TABLE 1 (CONI' 0) 



(Plus 4220 
& 4230, if 
appropriate) 


(Plus 4320 
A 4330, if 
appropriate) 


(Plus 4420 
& 4430, if 
appropriate) 


(Plus 4520 
A 4530, if 
appropriate) 










C. Data Read by Subroutine CINPUT. If I$TC?<Q 



^N0Z0N 


i NSTPZ 


IPUN1D 



lOENTIFICA TION 6112 73 


5 0 10 !^ 


IDENTIFICATION 6112 73 



1 t ». 


IDENTIFICATION 6E12.8 73| 


5 0 2 0 i80 

I ... t i _ . .A ^ 

APRgF(l,l) 

APRgF(l.2) 

.APRgF(2,2 ) 

APR0F(3,1) 

APR0F(3,2) 

5 0 3 0 




























TABLE 1 (CONT'D) 




(If NMR-12) 


(Plus 5220 A 
5230. if 
appropriate) 


(Plus 5320 fi 
5330, if 
appropriate) 


(Plus 5420 i 
5430, if 
appropriate) 














TABLE 1 (CONI' D) 


NUMBER 


61 i 


-t— ^ 


...i. 


* 


IDENTIFICATION fflg.a 


L'l 


- 1 — •- 






61 


1— — I ... -A., g 

-1 1 A A 


+ 


tDlNTtffCATWIt 6E12.8 


131 


-h 


■k, ^ I * 


IDENT I F ICATION 6E12.8 


_I3 

25 


- — 4 — ^ 


DESCRIPTION 


TMWm 


X=lUMMR- 




73 




. S.5J Q 1 


!80 


CSTR(I), 


I=1,NMR 


1 


■ 5,6,1 0 


80 


TVF(I) 


I^l.N TK 




73 


5710 j8P 


CP VF(I) 


I=KNTK 


(Plus 5520 A 
, 5530, if 
-1 appropriate) 


(Plus 5620 A 
5630, If 


! appropriate) 


! (Through 5740 
i f 

appropriate) 






IDENTIFICATION 6E12.8 73 


! (Through 58 « 
' if 


. ■ 5.8.1 0 80 


appropriate) 














TABLE 1 (CONI' D) 



NUMBER 

— 

OESCPIPTION ; 




TNBF 

111 



JSF .. _ 

- 

I 25 


^ A j j_ 1 1 . _k 1 t _ 

RH0N8F 

i37 


J. 1 1 ± j. . A. f * A 

RH0LF j 

49 



WTMLLF • 

1 

■ 

1 

WTMLVF i 


lOENUFlCATION 6E12.8 


80 i 

D 

■ 


"1 

TNB0 

i 

■ 

* ^ 1 . A __ 1 1_ _» A A A 

TB0 

25 



RH0NB0 1 

37 


_1_ 4 1 1 A . ^ .4. 

RH0L0 i 

49 


. ^ ^ ^ 1 ^ ^ . ! 

WTMLL0 1 

61 


1 

» ^ ^ A A ■ -A _ A _«_ _■! 


MTMLV0 ! 


1 

£J5 i2 d. 1.^ — : 

11 


_L _ . .L 1 .1 j 1 A A A 1 ^JCRITF_ _ _ 1 

.12 


. . 1 1 ^ ;TCR1T0 j 

III, 


. ■ . . I . . . . Idhvf J 

ill 


i 

DHV0 1 

49 


i - 1 t 1 1 1 1 t A 1 

^ ! 

j 

[ii 


1 

1 1 1 1 1 1 1 1 

A - - 



©ENTIFICATION 6E12.8 73 

XU- 

80 

r 


■■■ ■,-■■■- . ■— ,„■ , ■ , . 

1ST 

ll3 



NSSTI 

25 



NMSTI 

IS 

■ 

t 

ICRC 

49 


1 A 1 1 1 1 1 1 1 

IPRSST 

m 

■ 

1 

> 4 I 1 1 1 --i 


IPRHST 

IDENTIFICATION 6112 75 

.,.^.4 AO.. 

80 


51 





















TABLE 1 



IDENTIFICATtON 6E12.8 73 



‘ ‘ A — 


L.. „ 


IDENTIFICATION 5112 73 



i| ‘ ■ 


IDENTIFICATION C»12.8 731 


(CONT'D) 
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CRT^ 




1 

1 


I 



6. 5.8 0 


80 

















This is the end of the data deck if ILISPj^O. 

If IPUN3D=1, the data deck is completed with the punched 
card output from 3DC. 

REPRODUCIBILITY OF THE 
ORIUNAL PAGE IS POOR 



(Include cards 
only if ILISP=0 



S3 



TARLP, 1 (CONT’D'I 1-1 --4 
The following cards (7010 et seq.) are punched by subroutine STAPE if ILISP^O 
and IPUN|*0 in LISP input, as preparation for this input '»'her ILISP=0 



(Include cards 
only if ILISr=0 

and IPl'N?n=0'' 
i 


IXHUtj 


T 


54 
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TABL- 2. DER INPUT DATA DEFINITIONS AND UNITS 
(ALPRABETICAL ORDF«''l CARD NO. 


Fortran 
Variable Name 

Math 

Symbol 

Definition 

1 

„ • 1 
Units 1 

ALFA (I) 
lilSNEL 

APR0F(I,J) 

a 

The counterclockwise angle of element rotation 
about its z axis (see page 73) 

Combustion Chamber Geometry Specification 

^ degrees 

l^I^NAP 


APR0F(I,1) * Axial Distance 

in 



API^F(I,2) = Chamber Diameter 
If APR0F(1,1) ^ 0, APR0F(1,1) = throat 
radius ratio 

in 

A01(I) 

i<i2;lspec 


Cross-sectional area of orifice No. 1 of 
Type 6 or Type 9 element 

in^ 

AREAl(J) 

A 

Stream tube cross-sectional area at 

in^ 

iSjiNST 

z 

z - ZSTART 


ART0LD 

% 

Decimal tolerance* deviation of computed 
minimum stream tube area sum from nozzle 
throat area* multiple stream tube analysis 


BETA(I) 

l^I^NLSPEC 

3 

Counterclockwise angle of element rotation 
about its y-axis (see page 73) 

degrees 

CDDIA1,2(I) 

iSllNLSPEC 

S 

Discharge coefficient for orifice 1,2 of 
Ith element specification 

1 

CKP1,2 


Modified k* evaporation coefficient for 
propellaut l(fuel), 2(oxidizer) 
or Latent heat of vaporization (see page 78) 

1 2 
in /sec 

BTU/lbm 

CPVF,0(T) 

c 

Fuel, oxidizer vapor specific heat at 

BTU 

i^isntk 

pv 

constant pressure 

lbm-"R 

CRT0L 


Decimal tolerance, deviation of computed 
single stream tube throat contraction 
ratio from unity 


CST(I) 

c* 

Theoretical characteristic velocity at 

ft 

1<I<NCSTR 


mixture ratio 13MR(I) 

sec 

CSTR(I) 

1<I<NMR 

c* 

Same as CST 

ft 

sec 

DEAR-, 2 (I) 

kunlspec 

D 

Mass median droplet diameter for whichever 
propellant flows from orifice No. 1,2 of 
Ith e7^ment specification 

in 
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TABLE 2 (CONT’D 


Fortran 
Variable Name 

Math 

Symbol 

Definition 

^^nits 

DELTZG.L(I) 

I =. 6 or 10 

5 5 
G, L 

Axial location of pseudo-impingement point 
for gas (fuel), liquid (oxidizer)of 
Ith element specification (see page 72a) 

in 

{ 

t 

DENS 1,2 

0 

Density of propellant 1 (fuel), 2 (oxidizer) 
at injector temperature 

Ibra I 

^ ! 

DHVF, 0 

Ah 

V 

Fuel, oxidizer latent heat of vaporization 

BTU 

DIA1,2(I) 

1<I1NLSPEC 

D 

Diameter of injection orifice No. 1, 2 of 
Ith element specification 

in 

DNS ATI, 2 

P 

s 

Density of propellant 1 (fuel), 2 (oxidizer) 
at saturation temperature corresponding 
to the chamber pressure 

lb|n 

IP 

O^BAR(I) 

llIlND 


Spray size group diameter, D, divided by D ; 
with FRACUM(I) , gives spray distribution 

- 

DPINJ1,2 

Ap 

Injection pressure drop for propellant 1 
(fuel) , propellant ? (oxidizer) 

ibf 

in^ 

DRADM 

Ar 

Spacing between circumferertiai mesh lines 

in 

DTHETM 

A0 

Spacing between radial mesh lines 

degrees 

EPS 

ec 

Combustion chamber contraction ratio 


FRACUM(I) 

1<I<ND 

Mass fraction of spray assigned to Ith size 
group, whose diameter is given by D0DBAR(I' 

- 

GAME (I) 

llIlNLSPEC 

^E 

Stream impingement angle of Ith element 
specification 

degrees 

GAMFAN(I) 

1<I<NLSPEC 

Yf 

Fan cant angle of Ith element specification 
(Type 3 only) 

degrees 

GAMMA(I) 

1<I<NLSPEC 

Y 

Counterclockwise angle of element rotation 
about its x-axis (see page 7 3) 

degrees 

GASFL(J) 

1<J<NST 

w 

g 

Total stream tube gas flow at z^START 

Ibm/ sec 

GDIAD1(I,J) 

1<1<NGT 

1<J<NST 


Propellant spray size group droplet diameter 

I 

1 

in 

GVELD1(I,J) 
1< I<NGT 
1<J<NST 

“d“ 

Propellant spray size group velocity 

ft 

sec 
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TABLE 2 (CONT’D) 
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Fortran 
Variable Name 

Math 

Synbol 

Definition 

Units 

GWSPR(I,J) 

1<1<NGT 

liI:SNST 

. n 

"d 

Propellant spray size group flowrate 

Ibi 

sec 

ICRC 


Number of corrector cycles calculated at 
each Az interval (normally - 1) 

- 

IDBAR 


Integer^indicator of source of 5 data: 

D must be read in for all elements; 
=0, D calculated internally for NTYPEi.5 


ILISP 


Control integer; non-zero value causes 
subroutine LISP to be called 

- 

IPRMST, 


Number of Az intervals between multiple. 

-- 

IPRSST 


single stream tube printouts 


IPUN 

lPf1N5D 


Control integer; a non-zero value causes 
subroutine STAPE to punch-out stream 
tube initialization data- value of 1 

causes LISP to punch cards for 3DC 
Control integer; if equal to 1, STC reads 
punched cards from 3DC 


IPUNL,R 


Obsolete control integers; enter zeroes 

- 

IRCRT(I) 

l^ISNCRT 


Indices of circumferential (r=const) mesh 
lines along which CRT plots of mass 
flux are to be generated 


1ST 


Redundant control integer, set = 1 

1 

ISTC 


Control integer; non-zero value causes 
subroutine STC to be called 

- 

ITMC 


Obsolete control integer; leave blank 

1 

ITRANS 


Control integer, non-zero value causes 
subprogram TRANS to be called 

_ 

JSYM 


Integer controlling folding-in of flux 
from exterior radial mesh lines; 

« 0, no folding; = 1, mirror image 
folding; =2, repeating image folding 


KFCRT, 
K0CRT, 
KTCRT and 
KFFCRT 


CRT contour-plot controls, for fuel, 
oxidizer, total propellant and modified 
fuel fraction, respectively. A zero value 
suppresses the plot. A positive value 
denotes the number of contour lines. A 
negative value causes subroutine SCALE to 
select approximately that number of 
contours at rounded-off intervals. 

(Maximum values = 35) 
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TABLE 2 (CONT’D) 
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Fortran 
Variable Name 

Math 

Symbol 

Definition 

Units 

NRBAFL.R 


Number of circumferential mesh lines inter- 
sected by a radial baffle along the left, 
right boundary mesh line 

- 

NRML 


Number of circumferential (constant radius) 
mesh lines (S20)^ 

- 

NRWALL 


Number of circumferential mesh lines to 
the chamber wall 

- 

NSSTI 

! 

Maximum number of complete passes, marching 
from 2 * ZSTART to throat, in single i 

tube analysis 

■ 

NST 


Number of stream tubes(^40, but usual 17^10-20)^ 

- 

NSTPZ 


Number of stream tubes per zone^ually 4to 10 ) 

- 

NTHML 


Number of radial (constant 6)mesh lines(S20)^ 

- 

NTHL.R 


Index of the radial mesh line which forms 
the left, right boundary of chamber 
slice analyzed 

] 

NTK 

i 


Number of temperatures at which propellant 
vapor specific heats and film thermal 
conductivity are tabulated (^20) 


NTYPE(I) 

1<ISNLSPEC 


Element type index for Ith element 

specification (1<NTYPE^10) as numbered on 

page 70 . 

1 

1 ” 

NUG 


Index specifying uniform (NUG>0) or non- 
uniform (NUG = 0) gas velocities for 
gas/liquid elements 

1 

1 

PCI 


Static pressure at z * ZSTART 

1 psia 

PCTBL 


Percent (or decimal) of total propellant 
flow rate to be assigned to a chamber 
wall boundary layer stream tube 


PRPR1,2 


Propellant properties grouping (icj/e ) 
for propellant 1 (fuel) , 2 (oxidizer) 

f t-lbf 
sec 

PSTREC(I) 

liliNLSPEC 


Liquid injection post recess for a 
coaxial jet (Type 6) element 

in 

RADE(I) 

1<1<NEL 

'e 

Radial coordinate of the 1th element’s 
origin (e,g«, impingement point) 

in 


1. These limits are imposed hy the plotting routine PL0TC. 

2. The relationship NST = 1+N0ZON*NSTPZ holds if stream tubes are initialized 
via option (2), page 9. 
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TABLE 2 (CONT’D) 


Fortran 
Variable Name 

Hath 

Symbol 

Definition 

Units 

RH0G 


Combustion gas density at nominal chamber 
conditions 

^9 

RH0LF,0 


Liquid fuel, oxidizer density at saturation 
temperature corresponding to Pc 


RH0NBF.0 


Liquid fuel, oxidizer density at normal 
(latm.) boiling point 


RMN0M 


Nominal bulk mixture ratio 

- 

SA1,2(I) 

1£I<NLSPEC 


Spray coefficient a for orifice no. 1,2 
for Ith element specification 

- 

SB1,2(I), 


Ibid, coefficient b 

- 

sen, 2(1) 

to 

SC61,2(I) 

m^NLSPEC 


Ibid, coefficient 
to 

Ibid, coefficient 

- 

SMRG(J) 


Stream tube gas mixture ratio at 
z * ZSTART 

- 

SNN(J) 

1<J<NST 


Number of LISP mesh points combined to 
form a stream tube 

- 

SPEL(I) 

liliNLSPEC 

X 

Spacing between doublets of a like-Kioublet-* 
pair (Type 3) element 

in 

SPFAN(I) 

l^I^LSPEC 

Y 

1 Spacing between doublet’s spray fans of 
1 a like-doublet-pair (Type 3) element 

in 

SR(J) 

liJiNST 

i 

1 

1 

i 

1 Stream tube mean radius 

1 

in 

STH(J) 

liJ^ST 


Stream tube angular position 
! (set = 0 in subroutine SCRMBL) 

radians 

TBF,0 


Fuel oxidizer droplet saturation temperature 
at 

*R 

TC0NVF,O(I) 

k 

Thermal conductivity of vapor/gas film 

BTU 

1<I<NTK 


surrounding fuel, oxidizer droplets 

ft-sec-®R 

TCBITF,0 

T 

cr 

Fuel, oxidizer critical temperature 


TGAM(l) 

ISISNMR 

Y 

Combustion gas specific heat ratio 
(frozen) at TMR(I) 

- 

TGM(I) 

liUNCSTR 



Y 

Same as TGAM 
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TABLE 2 (CONT'D) 


Fortran 
Variable Name 

Math 

Symbol 

Definition 

Units 

TOi(I) 

T 

Combustion gas stagnation temperature at 

“R 

1<I<NCSTR 

o 

TMR(I) 


THETAE(I) 

l<r<NEL 

«E 

Angular coordinate of Ith element's origin 
(e.g.f impingement point) 

degree 

THETAR 


Angular coordinate of farthest-right radial 
mesh line 

degree 

TMR(I) 

l<IiNCSTR, 
or 1<I<NMR 

c 

Combustion gas mixture ratio array, defined' 
as the flowrate ratio of propellant 2 (oxi- 
dizer) to propellant 1 (fuel) 


TMU(I) 

1<I31CSTR 

U 

Combustion gas viscosity at TMR(I) 

ft=sec 

TMW(I) 

M 

Combustion gas molecular weight at 


lil^CSTR 

V 

TMR(I) 

mole 

or l^IiNMR 




TNBF.0 


Fuel, oxidizer normal boiling temperature 

’R 

TPR0P1.2 


Liquid fuel, oxidizer temperature used for 
initial vaporization calculation* Same 
as TBF,0 Is recommended but injection 
temperature may be used 

“R 

TT0(I) 

T 

Same as TG0 

"R 

li:XiaiMR 

o 



TVF,0(I) 

liI<NTK 


Temperature at which CPVF,0(I) and 
TC0NVF, 0(1) are tabulated 

"R 

ms (I) 
1<ISIMR 


Same as TMU 

Ibm 

ft-sec 

W0T1(1. 


Propellant 1 (gaseous fuel) flowrate for 


liliNLSPEC 


gas/liquid element specifications 
(Type 6 or 10) or flow for orifice 1 
of a Typft/9 element 


W0T2(I) 


Propellant 2 (liquid oxidizer) flowrate 


ISliNLS 


for gas/ liquid element specification 
(Type 6 or 10) 

sec 

WT‘tLLF,0 


Molecular weight of liquid fuel, oxidizer 

Ibn, 

mole 

WTMLVF,0 

1 

i 

Molecular weight of vapor fuel, oxidizer 

Ibm 

■ole 
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TABLE 2(C0NT’D) 


Fortran 
Variable Name 

Math 

Symbol 

W1F,0 


W2F.0 


NIFF 


W2FF 


WIT 


WlT 


ZE(I) 

Ifl^LSPEC 

*E 

Z0M 

^0 

Z0M2,3 


ZSTART 



Definition 


Lower limit of fuel, oxidizer flux contours 

Upper limit of fuel, oxidizer flux contours 

Lower limit of fuel fraction contours 

(= 0 .) 

Upper limit of fuel fraction contours 

(= 1 .) 

Lower limit of total flux contours 

Upper limit of total flux contours 

Axial coordinate of dementis origin (e.g*, 
impingement point) for Ith specification 

Axial coordinate of first plane in which 
spray mass fluxes are calculated 

Ibid, second, third planes 

Axial coordinate of stream tube initiali- 
zation plane 


1 . ZSTART may or may not be equal to the last Z^M analysed by LISP, 
as discussed on page 84* 
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In Table 1, the type of FORTRAN input (integer, floating point decimal, alpha- 
numeric) and the subdivision of each card’s first 72 columns into fields is 
indicated by its noted format. For decimal data, six 12- column fields are 
used; for integer data, either twelve six-column fields or six 12-coluinn fields 
are used. 

Standard FORTRAN input formats are used. Specifically used are: 

Comment Cards (A- format) 

Integer Variables Beginning with Letters I thru M 
(No decimal points, 12 (or 6) space field widths, 
last digit in last space of field, 6 (or 12 
consecutive values per card until READ statement 
is finished) 

Decimal Variables Beginning with Letters Other Than I thru M 
(Use decimal point or account for implied 
decimal location, one value every 12 spaces, 

6 consecutive values per card) 

The ’’Description” column in Table 1 gives the FORTRAN code names of input 
variables as they appear in the program listing. One value is to be entered 
for each coded variable unless it is subscripted. Array sizes for subscripted 
integer and decimal variables are indicated in this column below the variable 
names. (With A-format data, the subscript ranges are not indicated; 18 fields 
per comment card are implied.) For most of the data, all the values of one 
variable are read before proceeding to the next variable. (Exceptions are in 
the last group of stream tube initialization data which is part of the deck 
normally punched by STC’s subroutine STATE. The particular form used there 
permits use of very compact input/output statements and results in compact 
data decks, but is admittedly somewhat Inconvenient for manual generation of 
data decks.) 

3.3.1 Selection of LISP Input Data 

Organization of the LISP input data in Table 1 is, roughly, as six blocks of 
cards. The first block consists of two comment cards, in alphanumeric (A) 
format, which permit the user to document the case with such information as 


18A4 

6I12(or 1216) 


6E12.8 
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Injector name and/or drawing number, propellant combination, nominal chamber 
pressure and mixture ratio, date of the computer run, etc* The information 
provided on these cards is printed out with certain data printout headings 
and appears on computer plotted graphs* 

The second block of cards contains the description of the geometry of the 
system (number of elements, mesh lines to be set up, location of vail and 
baffles, if any, etc*), liquid propellant properties, and controls for the 
desired calculations and for CRT graphical output. 

The third block contains descriptions of common characteristics of groups of 
elements which make up the injector. Each of these groups is assigned an 
element specification index LSPEC which ranges from unity to NLSPEC. Data 
Included are the element type (doublet, triplet, etc.), orifice diameters 
and discharge coefficients, propellants issuing from the orifices, and 
geometric specifications, such as the axial coordinate of the impingement 
point and orifice alignment angles. 

The fourth block of data, which contains the information specific to indi- 
vidual elements, consists of NEL cards where NEL is the number of elements 
considered in the analysis. Each of the cards lists the specification group 
LSPEC which describes the given element, the radial and angular coordinates 
of the element's origin and an element orientation angle, ALFA, explained 
later* 

The fifth block of data, input if contour plots are called for, is concerned 
with the upper and lower limits of the contour lines. 

The sixth block of cards contains tables of the equilibrium combustion gas 

properties for the propellant combination. In the order of their read-in, 

they consist of the mixture ratio, the gas stagnation temperature (^R), 

the gas viscosity (lb /ft-sec), the ratio of specific heats, the mean 

m 

molecular weight (Ib^^iole), and the characteristic exhaust velocity (ft/sec). 
A minimum of two and a maximum of eighteen mixture ratios may bf put. 
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A functional, rather than a card^by-card approach is taken in the following 
discussions. Reasonably self-explanatory data are not commented upon. 

3. 3 .1.1 Geometric Considerations . Cards 2030 through 2060 contain a number 
of control integers and other data pertaining to the portion of chamber cross- 
section to be analyzed using an appropriate system of mesh points. LISP is 
structured to take advantage of injector symmetry, so that the data selection 
consists, briefly, of choosing a pie-shaped sector of the injector bounded by 
either planes of symmetry or baffles, such as shown in Figures 6 and 7, and C 
overlaying this sector with a grid of equally spaced radial and circumferential 
mesh lines. In both Figures 6 and 7, the area bounded by AOC becomes a satis- 
factory segment for analysis. LISP is dimensioned for up to 400 mesh points 
which permit the choice of combinations of as many as 20 x 20 or 25 x 16 mesh 
lines if desirable. Similarly, LISP is also dimensioned to describe up to 
60 injector elements. 

It is often appropriate in LISP to employ additional mesh lines and injector 
elements outside the chamber slice under analysis. The spray flux within the 
slice to be analyzed includes spray from nearby injector elements when the 
sides of the slice are planes of symmetry; spray from interior elements will 
collect on the baffles when they constitute the side walls; and some spray 
will Impinge on the walls of the thrust chamber ahead of the plane being 
analyzed and /un down the walls to add to the spray directly impinging on 
wall mesh points in the plane. To treat these conditions, a number of options 
are included in LISP. The options can be described by reference to the sample 
mesh line system depicted in Fig. 8, which represents an overlay of a mesh 
line system on the injector segment of Fig. 7. In this illustration; 

1. NRML* is the number of circumferential mesh lines. 

2. NRWALL is the number of circumferential mesh lines to the chamber wall. 
When (as in Fig. 8) NRWALL is less than NRML, the spray mass to the mesh 
points beyond NRWALL is folded into the wall at the corresponding 
6-location. 


^^e nomencTat'ure for F(?RTRAN variables is given in Table 2 
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Figure 7. Selection of Segment of Injector for Analysis of 

Spray and combustion Gas Flow Fields - Segment Bounded 
by Combination of Planes of Symmetry and Baffles 
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Figure 8. Overlay of bystem of Mesh Lines Upon injector 
Slice to Permit Calculation of Spray Flux 
Distribution by LISP 


2=NTHR 

1 = 1 

NRML 

J KJR3AFR 



3. NTHML is the total number of radial mesh lines. The radial mesh lines are 
indexed J*1 to J-NTHML counterclockvise (i.e., from right to left). 

4. NTHR and NTHL are tne radial mesh lines which define right and left 
boundaries of the chamber slice to be analyzed. It should be noted that 
additional radial 15nes are defined (in this case) to the right and left 
of NTHR and NTHL, respectively. 

5. Because a baffle corresponds to the right boundary of the slice, a code 
number NRBAFR is input to LISP, which corresponds to the nuad>er of cir-* 
cumferential mesh lines (starting from the center) over which the baffle 
extends. (Since the baffle extends to the wall in Fig. 8 ^ NRBAFR is equal 
to NRUALL in this example.) The spray flux to mesh points to the right 

of NTHR wi31 be folded into NTHR at the corresponding radial locations. 
Similar results can be produced for baffles on the left boundary of a 
slice by the code number NRBAFL. The baffles must begin at the center of 
the chamber but they need not extend to the wall, i.e., NRBAFR and NRBAFL 
may be less ban NRWALL. 

6. The best a*. .i erred way to account for the fact that the spray inside a 
slice may include contributions from nearby elements outside the slice is 
simply to include these ^elements in the input specifications. For example, 
in Fig. specifications for til elements within the dotted lines would 
be made in the analysis of the slice AOC. 

7. The procedure described in (C) above is the reconmended method 
accounting for spray travel across slice boundaries defined by planes of 
symmetry; however, alternate methods are provided in LISP by the F0RTRAN 
code variable JSYM^ When JSYM*0, the procedure in (6) above is employed. 
When JSYM is defined as either 1 or 2, only injector elements within the 
area encompassei by radial mesh lines NTHR and NTHL are input as data, 
but the spray from the elements to mesh points outside the area are folded 
back into the corresponding mesh polices inside the area. When JSYNhl, 
this folding is done on the basis that the elements to the right and 

left of both the boundaries NTHR and NTHL are mirror images of each 
other. In this case, the spray flux to mesh line (NTHR-1) is folded 
back into meah line (NTHR-f 1) ; similar relations hold around boundary 
NTHL. Ii the elements in and around slice AOC are related as repeating 
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nets, JSYMa2 is defined and the folding of flux from exterior to interior 
mesh lines is done according to (ifFEBUl) folds into (NTHL-l), etc. 


In analyses of aereral injector designs at Rocketdyne« it has been experienced 
that choosing JSYM-0 offers fearer "trapa** in setting up a probleai, particu- 
larly where the geometric distribution of injector elements may require 
assusiptions of effective planes of sysmwtry which are not rigorously so. 
Although many plcuiea of syamMtry may be defined for an injector, analysis 
of the gas flow within a given chamber slice will give meaningful results 
only if the injection elements are specified such that the correct geo- 
metric proportions of total propellant flows enter the slice; otherwise, 
erroneous mass flux levels and chamber pressures will be calculated. For 
JSIMbO, sillily defining enough elements outside the chamber slice takes care 
sf this. For JSYM>lor 2, elements that fall on (or very near) the 
boundaries should be examined carefully and their hole sizes or flow rates 
may have to be adjusted artificially to ensure proper flows. 

3 .3. 1.2 Injection Element ConSi. derations . The current version of LISP permits 
consideration of the following t«i t 3 rpes of injector elements: 

1) Ikilike doublet 

2) Like-doublet 

3) Like-doublet-pair 

4) Triplet 

5) Four-on-one 

6) Gas/liquid coaxial jet or concentric tube 

7) Showerhead* 

3) Special callout by general spray flux equation 
9) Gaseous transpiratioi cooled injector face 
lO) Gas/liquid triplet 


*Distribution coefficimts have not yet been evaluated for the showerhead 

element, which must currently be treated as a Type 8 element. 
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Element Specif icat ions ♦ Card 2030 carried a value of NLSPEC^ the number of 
element specifications which may range from 1 to 10. If all elements are alike» 
differing only in their position and orientation on the injector face» a single 
specification is needed (NLSPEC*!)* Any element or group of elements that differs 
from the others in element type, element design (e.g. , orifice sizes, discharge 
coefficients, impingement angle and distance), propellant assignment to the 
orifices or element angular alignment (other than a simple rotation about its 
z-axis) requires a separate element specification. 

Element spe Ification data ore entered in cards 2x10 through 2x95 where the 
X denotes eac’ specification number* Mass flux distributions calculated by 
means of Eq* (1) and (2) make use of correlation coefficients built into function 
subprograms in terms of such parameters as the element types, orifice sizes, 
stream impingement angles and momentum ratios, etc*, for Type 1 through 6 
elements* The correlation coefficients, must be supplied by the user for 
Types 8 and 10 elements* For the first five element types listed above, the 
element flow rates of each propellant are computed within the LISP subprogram 
EFLOW in terms of the element orifice sizes, discharge coefficients, and the 
injection Ap*s of the separate propellants* For type 6 and 10 gas/liquid 
elements, propellant flowrates are supplied as input data. Mean drop sizes 
for propellant sprays may be supplied directly as input data for all elemental 
types or, in the case of the first five element types, may alternatively be 
calculated within the LISP subprogram DSIZE. 

If an element or group of elements is defined as being type 8, the mean droplet 
si?3 is read in as input data and the total flow from the element is calculated 
as if the element were a doublet with equivalent orifice diameters. The mass 
flux distribution is calculated by means of Eq. (2) using spray distribution 
correlation ccefficients supplied directly as input data rather than being 
supplied by subprograms %d.thin LISP. This feature permits a designer to have 
a cold-flow characterization made of the single element (or elements) to be 
incorporated into a prospective injector and then to employ the correlated 
spray coefficients from the cold flow experiment in LISP. Such a procedure 
is useful .ot only for the situation where spray coefficients have not been 
determined previously for th? intended elements, but it al9o permits the designer 
to acco r Tor factors such as short L/D orifices and m-s^iifold cross-flows in his 
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single element cold-flow experiments. A Type 8 callout may also be used to lump 
together into a single effective source a number of elements near the center of 
a large injector. 

In the ease of a transpiration-cooled injector face, the Type 9 element flow 
rate is read-in and is divided by the cross-sectional area of the cliasber 
slice to define a uniform mass flux at all r, 9 mesh points. 

Single Flement Geometry . LISP converts the (r,9,z) coordinate system of 
the thrust chamber into individual rectangular (x,y,z) coordinate systems for 
each eleoient based upon particular origins of the elements. Element coordinate 
systems are illustrated in Fig. 9. For the liquid/liquid doublet and triplet 
elesmits, the x axis corresponds to the long axis of the spray fan formed by the 
impinging streams while the y axis lies in the plane of the impinging streasis. 

For the 4-on-l and like doublet pair elesients, the geometries shown in Fig. 9a 
were made as coosistent with the doublet and triplet as was conveniently possible. 
The z-axis origins for liquid/liquid eleswnts are at the element impingemnt 
points and, so far as input data are concerned, are always directed along the 
bisector of the stream impingement autgle. (For unlike doublets, the z-axis is 
changed intemall' by LISP to lie along the resultant spray momentun rector.) 

The definitions cf x and y axes for the type 8 element are equivalent to those 
specified for the liquid-liquid unlike doublet but the z-axis is not altered to 
follow the momentum vector. From synmetry, the definition of x and y axes for 
showerhead elements is inmaterial unless they are canted with respect to the 
chamber axis. 

The coordinates for gas/lir id triplet and coaxial jet elements are illustrated 
in Fig. 9b. The triplet definitions parallel thotu- for the liquid/liquid 
triplet, with an origin at the geometric impingement point. Coaxial jet eleisents 
are axisymmetric so that, like the showerhead, definition of x and y axes 
amtters only if an element is canted with respect to the chamber axis. The 
origin for the coaxial jer lies on the eleisent centerline and at the discharge 
end of the liquid (oxidizer) injection tube (post). If there is a "post-recess," 
the variable will be negative. 
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Figure 9a. Liquid/Liquid Injection Element Designations 
and Coordinate Systems 
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Figure 9b. Gas/Liquid Injection Element Designations 
and Coordinate Systems 
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LISP considers both the liquid and gas injected with a gas/liquid element to be 
distributed as if they flowed as rays from point sources* Rather than lying at 
the element’s coordinate origin^ as with the liquid/liquid elements, the point 
sources are located at distances of and upstream of the origin for the 
liquid and gas, respectively* These parameters are derived during correlation 
of elemental cold-flow distribution data and values must be read-in for each 
gas/liquid element. 

The orientation of the x, y, z coordinate system of an individual element to the 
chamber coordinate system is defined in terms of three rotation angles whose 
respective Fortran names are ALFA, BETA, and GAMMA* If; 

1* The elraent is oriented on the injector such that its y axis (as defined 
by Fig. 9 ) is coincident with the chamber radius through the element 
impingement point, with positive y pointing radially outward, and 
2. The z axis of the individual element (as defined in Fig. 9) is parallel 
to the thrust chamber axis , 

then the angles ALFA, BETA, and GAMMA have zero values. If the element is not 
oriented in this "basic** or reference alignment, then 

!• ALFA is the counterclockwise angle of rotation around the basic z axis 

of the element *s y axis from its original alignment with the chamber radius, 

2, BETA represents counterclockwise rotation around the y axis (in its new 
position after the first rotation); finally, 

3, GAMMA represents counterclockwise r*^ cation around the x axis (in its 
transformed position after the first two rotations). 

For each rotation, ’’counterclockwise** implies the rotation direction that would 
be seen by an observer looking along the positive axis toward the origin of the 
element. In most applications, no more than one rotation will be applied to 
the geometry cf any single element. 

Each element, except for Types 7 and 9, is considered to consist of two 
equivalent orifices designated as 1 and 2. The physical orifices which cor- 
respond to these number designations are listed in Table 3 for all element 
types. 
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The LISP program also labels the injected propellant as either 1 or 2, The 
fuel must be chosen as Propellant 1 and the oxidizer ..s Propellant 2 because 
of the combustion gas properties vs mixture ratio tables and because they are 
expected to bear these designations when data are transferred from LISP to 
STC« Also 9 for a gas/liquid propellant confcination, the gaseous propellant 
is always assumed to be the fuel. 

Propellant Flowrates . Except for injector elements defined as T 3 rpes 6, 9 
or 10, the weight flow Ib^/sec) of each propellant from each orifice of the 
element is calculated in LISP by means of the standard orifice equation. The 

3 

LISP user must input the injection density (lb /ft ) and the injection pressure 

m 

drop (psi) for each propellant as well as the diameter and discharge coefficient 
for each orifice of the element. If the injector element is specified as 
Type 8, then the user must specify the diameter of the circle of equivalent area 
as the orifice diameter because LISP treats the Type 8 element as an unlike 
doublet. 

The T}rpe 6 gas/liquid coaxial jet element flowrates are defined as follows* 

Liquid oxidizer flow is calculated from element specification data on post 
diameter, DIA2, discharge coefficient, CDDIA2, and injection pressure drop, 
DPINJ2. The gaseous fuel flowrate. W0T1, is read-*in along with the actual 
cross-sectional area, A^l, of its atmular injection passage. 
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Table 3* Physical Orifices Corresponding to Equivalent 
Orifice Designations One and Two for the 
Element Typ^s Defined in LISP Program 


Element Type Equivalent Orifice 1 


(Unlike Doublet) 


2 

(Like Doublet) 


3 

(Llke-Doub le t-Pai r ) 

4 

(Triplet) 

5 

(4-on-l) 

6 

(Gas /Liquid Coaxial) 

7 

(Shower head) 

8 

(Special Callout) 

9 

(Transpiration Face) 


The orifice lying on the 
negative y axis of Fig* 9 
which injects toward the 
positive y axis- It is the 
fuel orifice of Fig. 6 
(where the r coordinate 
corresponds to the y axis) . 


Same as for the Type 1 
element; however, the 1 and 
2 designations are not 
important for the Type 2 
element. 


The pair of ports to the 
left of the y-axis of Fig. 9. 


The two outer impinging 
orifices of Fig. 9. 


The four outer impinging 
orifices of Fig. 9 


The annual orifice through 
which gas flows* 


The single port. 

Treated as unlike doublet. 

The entire cross-sectional 
area. 


10 

(Gas/Liquid Triplet) Same as Type 4 


Equivalent Orifice 2 


The orifice lying on the 
positive y axis of Fig. 9 
which injects toward the 
negative y axis. It is the 
oxidizer orifice of Fig. 6. 


Same as for the Type 1 
element. 


The pair of ports to the 
right of the y-axis of Fig. 9 

The inner orifice of Fig. 9. 

The inner impinging orifice 
of Fig. 9. 

The inner post orifice 
through v;hich liquid flows. 

Not defined. 

Treated as unlike doublet. 

Not defined. 
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Diacharg* coafficianta of 0.84 are recommended for liquids flowing through short 
L/D orificaa under motor firing conditions and 0,80 under cold flov^ condition. 
However^ the user is free to use whatever coefficients he deems appropriate. 

He may^ for example^ use different coefficients within an element to account 
for effect of manifold cross flows in the injector, or use different coeffi- 
cients from element to element to approximate injector maldistribution effects. 

Spray Drop Sizes * If the injector element is specified as being Type 1 
through Type LISP will calculate a mass median drop diameter for the propel- 
lant of each orifice of the element. These calculations are based upon the 
correlations of Dickerson et. al. (Ref . l2) derived from hot wax experiments. 
Constants in the correlations have been modified to give characteristic 
diameters which make calculated c* efficiencies compatible with measured results 
for three injectors tested, analyzed and reported in Ref. 5. Alternatively, 
by assigning a value greater than zero to the flag IDBAR, the LISP user may 
assign his own estimation of drop diameter to the flow from each orifice of a 
given element* For elements defined as Types 6 through 10, the user always 
supplies his estimation of a characteristic drop size for each orifice of each 
element* The appropriate mean droplet diameter is the mass median diameter. 

With elements of Types 2 and 3, the mean drop diameters are based on the 
empirical correlation of Falk, etc. (Ref. 13), modified to make c* efficiencies 
calculated by the STC computer program correlate with experimental data reported 
in Ref. 13* Additional input data are required for these elements: combustion 

3 

gas density (Ib^/ft ) at injection-end chamber pressure, chamber contraction 
ratio and, for each propellant, the liquid physical properties grouping ya/p 

3 

entered as the parameters PRPRl and 2 on card 2080. Here p is density (Ib^/ft ), 

a is surface tension (Ib^/ft) and y is viscosity (lb /ft-sec), all evaluated at 

X m 

injection conditions* 

Propellant and Combustion Gas Properties * The densities of both propellants 
at their injection temperatures and at their saturation temperatures corresponding 
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to the expected chamber pressure are required* (Cards 2070 and 2080) • The 
former are used to calculate injection flow rates and momenta; the change in 
density to the latter are used in defining the output droplet diameters* 

Propellant temperatures (Card 2090) are used in estimating an appropriate value 
of the evaporation coefficient, Normally, these should be supplied as the 

saturation temperatures corresponding to the expected chamber pressure. 

The required combustion gas properties are supplied as tables which relate the 
given gas properties to mixture ratio at assumed chemical equilibrium conditions. 
The tabular property data supplied are lists of the stagnation temperature,** 
viscosity, molecular weight, specific heat ratio (y) » and characteristic exhaust 
velocity (c*) corresponding respectively to a list of mixture ratios. These 
data are given on Cards 4110 through 4630. These combustion gas property tables 
have been set up so that essentially the same set o^ punch cards may be used 
again as STC input data. 

The number of entries in the gas property tab?es are defined by the control 
variable NCSTF. A maximum of eighteen mixture ratios may be used in the tables 
which should range from essentially all fuel vapors (zero or near zero mixture 
ratio) to essentially all oxidizer vapors (a mixture ratio of at least 100) to 
prevent possible extrapolation errors when local regions of extreme propellant 
maldistribution car encountered in LISP analyses. If the LISP user is inter- 
ested in calculating only the spray mass distribution for a given injector 
configuration, then only two entries (one punch card) are required for each of 
the six gas property tables listed above (NCSTR*2). 


*In the case of a gas/liquid propellant combination, the density of the gaseous 
propellant is supplied at its injection temperature and the expected chamber 
pressure. An error in the estimation of chamber pressure v/ill effect only 
initial vaporization and dropsize and will not effect the definition of 
propellant flow rate. 

**The required properties may be obtained, for example, from runs of the 
JANNAF 0DE (one-dimensional kinetic) computer program (Ref. 14) or an 
equivalent program (e^g., Ref. 15) 
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Partial propellant evaporation between its injection and the LISP collection 
plane is approximated by Eq. (3). The user needs to concern himself with a 
value of Cj^, for each propellant. One method is for the user to input values 
of Cj^, for the variables CKPl and 2 on Card 2070; the obvious implication is 
that he knows appropriate values. If that is not the case, another method has 
been built into LISP; values of are calculated internally if the user 
inputs values of latent heat cf vaporization for CKPl and 2. The program 
distinguishes between these methods by examining the magnitudes of the input 
values: if both are larger than unity, the values are interpreted as heats of 
vaporization. Conversely, if either is less than unity both are retained as the 
modified evaporation coefficients. 

If CKPl and 2 are heats of vaporization, are calculated in subroutine 

CCKP. An empirical scale factor has been applied to the theoretical values to 
force agreement between LISP's calculated percentages burned at Zq» 2 and those 
observed in an NT0/N2H^-UDMH (50-50) engine experiment (Ref. 5). 

_A 

Numerical values of C. , will ordinarily range from approximately 2 x 10 to 
-3 2 

2 X 10 in /sec. Typically, they are on the order of 1/5 to 1/3 the values 
computed later, in STC, for for the same propellants. This is appropriate 
because Cj^, is being applied, in LISP, to a spray that is incompletely atomized 
and whose droplets are not yet heated to their saturation temperatures over much 
of the travel distance considered. 

CRT-Plottine Options . Input data concerning computer generation of data 
plots are scattered throughout the LISP input data. A non-zero value of NCRT 
on Card 2030 will result in a cross-sectional plot of the LISP mesh system, 
with all element origin locations denoted. If NCRT is larger than zero, tha«. 
number of mass flux profile plots will be produced, each showing the fuel and 
oxidizer fluxes vs angular position, 9, at a fixed value of chamber radius. The 
radial coordinates for these plots are prescribed through the variable array 
IRCRT on Card 2050 (et seq.). Each value of ICRT denotes ..ije number of a cir- 
cumferential mesh line, counting from the center of the chamber out. 

Contour plotting is controlled by the values of the variables KFCRT through 
KFFCRT on Card 2040. A non-zero value for any of these variables results in one 
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contour plot. A positive value denotes the number of contour lines which are 
to be plotted; value3 between 10 and 15 ordinarily produce pleasing graphs. 

A negative value denotes only approximately the number of contour lines; the 
actual number is determined internally by selecting a rounded-off contour 
interval. Use of a negative value is recommended. 

The maximum and minimum ranger of contour plotted variables can be controlled 
by input data on Cards ^010 and 4020. (These cariJ- ' e needed if any of 
KFCRT through KFFCRT are non-zero.) If limit values are input as zero, 
contour scaling is done internally. Ordinarily, one would not know, a priori, 
appropriate ranges for mass fluxes, so zeros would be entered for WIF, 0 and 
T and for W2F, 0 and T. On a later run, it might be desired to examine some 
contour interval in more detail by excluding some higher or lower flux values. 
Ordinarily, the modified fuel fraction is restricted to the range W1FF=0, 
W2FF«1. 

Subscripts . The variables which appear in the LISP calcu*’ itions are sub- 
scripted according to either injector element index (1-60), combustion zone 
mesh point index (1-400), propellant index (1-2), orifice index within a 
specific element (1-2), or element spcciKcation set (1-10). Of these indices, 
the element index, the mesh point index and the element specification indices 
become the regular Fortran integer indices I, J, etc., for use in D0 loops '»nd 
in READ and WRITE statements. The remaining subscripts are incorporate’ 
directly into the Fortran variable namcc.. For example, the weight fh c 
Propellant 1 from Orifice 2. of injector element 3 would appear in the Fortran 
code as W0T21(3) and spray flu a of Propellant 1 from all elements to mesh 
point 22 would appear as the Fortran variable STW1(22). Interpretation of 
the program listing in Volume 2 and the input variable designations included 
in Tables 1 and 2 should be made in this context. 

Multiple Collection Plane Analyses . By specifying non-zero values for 
Z0M2 (and Z0M3) , spray mass dictributions will be calculated by LISP at a 
second (and third) collection-plane location. Data transferred via scra». ’h 
data unit 2 to STC are from the last collection plane analyzed. 
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3.3»2 Selection of STC Input Data 

Input data for the STC subprogram block are read by subroutine CINPUT. Card 
sequence numbers begin with 3010 » which provides values for some stream tube 
initialization control parameters. If ILISP*0, stream tube initial data are 
read-in via punched card so the parameters NSTPZ and NUG have no meaning. 

In that case 9 the value of N0Z0N determines whether the stream tubes are 
analyzed assuming axisymmetry (N0Z0N>O) or not (N0Z0N*O) . 

If ILISPf^O, stream tubes are to be initialized from LlSP-generaiied mesh point 
flow data. Then N0Z0N-O is forbidden and will cause job termination. N0Z0N~1 
results in all LISP mesh points ou each circle of mesh re, s those 

at equal chamber radii) being combined into one stream tube (option 1, page 8) 
and NST*^NRWALL. If N0Z0N>1, its value denotes the number of circular or 
annular chamber zones (option 2, page 9) of mesh points a NS'^PZ gives the 
number of stream tubes per zone. Then, recognizing that the wall mesh points, 
at least, arr combined into a wall-bounding stream tube, the total number of 
stream tubes is NST»N0Z0N*NSTPZ + 1. 

With ILISP^O, a non-zero value of NUG causes all of LISP’s gasified propellant 
flow to be distributed uniformly among the mesh points (i.e., constant gas 
flux) • 

The option to analyze stream tubes by suppressing data concerning geometric 
mesh point locations and forming stream tubes entirely acc rding to a mixture 
ratio ranking is given by ILISP = 0 and a 0. Because LISP will be 

bypassed in this case, the portio*' of the data deck which initializes stream 
tube data must be provided (e.g*, from a previous run iu which LISP was run) 

If a non-axisyrametric multiple stream tube analysis is desired, zero values 
should be input IL'SP and ITRANS. Also, it should be recognized that only one 
multiple stream tube analysis will be performed in this case, i.e., there will 
be no iterative analysis to converge on a solution satisfying a throat boundary 
condition. This does not prevent, however, the punchout of throat-plane data 
in NAMELIST form for subsequent input to improved TDK computer program. 


80 


Card No> 5020> The variable NP deaign««tes the number of z-planes between 
Zq and the nozzle throat (inclusive) at which program calculations are to be 
made* The preferred method is to select a desired magnitude for Lz and 
calculate NP by: 


NP 


Az 


+ 1 


It is recommeiided that values of Az be selected in the range 0*02 <Aa^ 0*10-in* 
such that there are some convenient integral nuiri>er per inch for improving 
v'rintout readability* The other variables on this card are self-explanatory* 


Combustion Chamber Geometry (Card No* 5030* et seq*) * The geometry of the 
cony>ustor is described through the doubly subscripted array APR0F(J»L)* This 
array is entered in matched pairs (L»l,2) for each value of J* The values 
APR0F (J»l) denote axial distances (in) from the injector and the values 
APR0F (J,2) denote the corresponding chamber diameters (in) at these positions* 

It is required that distance increase with increasing J and that the array 
progress from the injector plane to the nozzle throat* That i$» APR0F (1»1) 
is assumed to be zero and APR0F (NAP^l) is the injector to throat distance* 

The intermediate values with KJ<NAP are used to describe the wall profile 
of an axisymnetric chamber* For z-planes lying between any APR0F (J,l) and 
APF0F (Jfl^l)^ linear interpolation on diameter is used* 

A recommended option is provided for cases employing the common design of a 
conical nozzle convergent section which is tangent to a radius of curvature 
through the throat section* This option is invoked by entering the value 
of the throat radius ratio (wall curvature/ throat opening) as APR0F (1,1)* Later^ 
in subroutine AVAR, it is reset equal to zero for subsequent interpretation as 
the injector face plane. Then, the upstream end of the nozzle cone is specified 
by the next to the last point, APR0F (NAP-1, L), and the nozzle throat by last 
point APR0F (NAP,L)* The nozzle profile between those two points is computed 
within subroutine AVAR* 


There are other options available in subroutine AVAR which are not expected to 
be useful for analyzing axisymmetric cylindrical chambers so their influences 
on APR0F (J,L) input are not discussed here. 
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Multiple stream tube STC analysis is continued past the throat for several 
z-planes« With the throat plane denoted by Z(NP), mirror-image symmetry is 
assumed such that areas at Z(NP+1) and Z(NP-l) are equal, etc, 

CoBfeustion Gas Properties (Cards No, 5130, 5210-5630) . Combustion gas 
properties, tabulated as functions of gas mixture ratio, are obtained from 
prior, peripheral computation using a thermodynamic equilibrium performance 
computer program, Rocketdyne’s free-energyninimization program (Ref, 15) 
has been used here, but any comparable program is adequate. Data entered in 
this table are properties for equilibrium combustion products at stagnation 
conditions corresponding to the mean expected chamber pressure. The 
parameters are assumed not to vary appreciably with pressure, but are taken 
to be functions only of mixture ratio. 

Propellant Vapor Properties (Cards No, 5710-6230) , Tables of fuel and 
oxidizer vapor specific heat and thermal conductivity as functions of tempera- 
ture are required, spanning the range of temperatures across the vapor /com- 
bustion gas films around spray droplets. At the lower temperatures, specific 
heats at constant pressure may be conveniently obtained directly from propel- 
lant enthalpy tables or charts* At higher temperatures, dissociation is 
important and it is ap cpriate to blend the low temperature, undissociated 
data into equilibrium dissociation data. 

The thermal conductivity needed is not simply that of the vapor, but that of 
the combustion gas-vapor mixture between a droplet’s surface and a surrounding 
flame-front. Again, a blending between undissociated propellant, dissociated 
propellant and propellant-rich combustion gases is appropriate at low tempera- 
tures with a gradual shift to the conductivity of the combustion gases alone 
at high temperatures. 

These properties are used In subroutine KPRIME to calculate propellant droplet 
evaporation coefficier.s, £q, (23). More detailed discussions of these data 
concerning selection and the following miscellaneous properties) appear in 
Reference 16 together with some example data. 
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Miscellaneous Propellant Properties (Cards No, 6510 to 6330 ) « The various 
temperatures, densities, and molecular weights are explained adequately in 
Table 2, as are the liquid molecular weights* The vapor molecular weights 
will normally be lower than those of their liquids; values corresponding to 
equilibrium dissociation at the mean film temperature are recommended here. 

Such dissociation is considered, Iiowever, to take place far enough from the 
droplet surface that it has little effect on the thermal gradients there. 

Hence, it is no longer deemed appropriate to add a heat of dissociation to 
the heat of vaporization. Similarly, approximating a droplet warm-up effect 
by adding a liquid sensible heat to the heat of vaporization is no longer 
recommended, i.e., the value entered for latent heat of vaporization should 
be just that. 

The heats of vaporization are used in calculating evaporation coefficients, 

Eq. (23). The liquid saturation densities are used in droplet size and 
behavior calculations and the rest of these properties are used in obtaining 
mean film viscosities in subroutine EVAPS. 

Program Control Variables (Cards No. 6540 and 6530) . An averaged, single 
stream tube analysis of the multiple stream tube input is forced by internal 
assignment of IST-1, regardless of the value read-in. The main purpose of 
this is to ensure that the initial-plane chamber pressure (PCI) and total 
propellant flowrates are compatible with the nozzle throat boundary condition 
of maximum gas flux. The boundary condition Is approximated by requiring 
that the single stream tube’s gas velocity at the throat be eoual to calcu- 
lated sound speed for the combustion gases at the calculated throat mixture 
ratio. A required flow area is computed from a gas continuity equation; the 
ratio of this area divided by geometric throat area may differ from unity 
by CRT0L or less; otherwise, the value of the initial pressure is adjusted 
and the entire single stream tube analysis is repeated. Convergence of the 
area ratio to (1.00 t CRT0L) is required within NSSTI iterations before 
proceeding to multiple stream tube analysis. Values of CRT0L = 0.005 and 
NSSTI = 5 are recommended. 
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Single stream tube data arc vritten out after each of the first four Az 
increments. Thereafter, the data are written out at intervals of IPRSST Az's 
until the neighborhood of the throat is reached. Data are then vritten out 
at every increment. Priit-out control in multiple stream tube analysis is 
identical except that the intermediate print interval is IPBMST. 

Improved computational accuracy may result from repeating the sequential spray 
phase/gas phase solutions using updated values of flow parameters, in a 
"corrector" cycle. Repetitive corrector cycles may be obtained by means of 
the variable ICBC which denotes the number of corrector cycle calculations 
to be made in each Az increment. The recommended value is ICRC = 1. 

The FORTRAN variable ART0LD denotes the decimal tolerance 6* discussed on 

A ^ 

page 18. A value of AHT^LD = 0.0025 has been used successfully in checkout 
analyses. 

Other Data (Cards 6560-6610) . The ZSTABT plane is used as the initial 
plane for vaporization calculations in STC* Its location is not necessarily 
the same as Z^, the LISP spray flux collection plane. The pattern calculated 

for the last value of Z0M used in LISP is carried over into STC. Then, if the 

value of ZSTART is in the range 0<ZSTART<Z^M, it may differ from Z0M. However, 
if ZSTART a 0 or ZSTART>Z^M is entered, it is reassigned as being equal to Z(Im. 

The variable PCTBL indicates the minimum percent (or, equivalently, decimal 
fraction) of the total propellant flow rate to be assigned to the stream tube 
adjacent to the chamber wall (if ILISP/o). A zero value is recommended unless 
the STC computed data are going to be used in a wall boundary layer analysis. 

The two arrays FRACUM and D^DBAR allow the user to control the distributions 

of spray droplet sizes about the values of B obtained from LISP, If these 
variables are not provided, the empirical distribution observed for like- 
doublet elements (Ref. 12) is used by default: 


Cumulative 
Mass Fraction 

0.05 

0.15 

0.25 

0.35 

0.45 

0.55 

0.65 

0.75 



D/B 

0.375 

0.57 

0.69 

0.8 

0.93 

1.07 

1.22 

1.39 

1.6 

2.1 
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Stream Tube Initialization Data (Card No> 7010> et seq.) Stream tube 
gases are completely defined by input of pressure^ flow area* flow rate and 
mixture ratio at a starting plane. All other desired parameters are calculated 
from these using the combu- tion gas properties tables. The spray flow rate, 
velocity, and droplet diameter for each spray group define the liquid input. 

Of the NGT spray groups, th^ first NGF are fuel spray and the remainder are 
oxidizer spray. 

Card 7010 carries the area, gas flow rate and mixture ratio for the first 
stream tube. Additionally, it has fields for the number of mesh points that 
went into the stream tube and the mean r and 6 coordinates of the stream tube. 
These last three variables are holdovers from an earlier program: they are 
written out but are not used in STC analyses. It is thus immaterial vhat 
their values are* 

3.4 PROGRAM OPERATION 

In addition to understanding program input and output, a program user needs 
to know operational information, such as required core size, typical execution 
times, reasonable line count and CRT limits and required auxiliary data 
units. Experience in these regards on Rocketdyne*s IBM System 360, Mod 50/65 
will be related briefly to provide preliminary indications, at least. 

With the overlay structure shown in Fig. 5, the program length is 42,450 
words. Total storage used during execution was about 47,200 words. 

It is difficult to be quantitative about execution times for several reasons. 
Foremost among these is the 360 *s ability to process several jobs simul- 
taneously so that **time** becomes rather nebulous. Neither clock time nor 
CPU time are valid indicators of what it costs to run cases. For that 
reason, figures used here are Rocketdyne ^billing units”, which are calcu- 
lated by a complicated formula but may be thought of roughly as "minutes of 
execution time”. 

LISP runs typically cost about 1-1/2 to 3 BU (billing units) for analyzing 
one collection plane. Higher costs go along with many injection elements 
and many mesh points. 
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STC runs typically cost 3 to 8 BU. Initialization and single stream tube 
analysis probably accounts for about 1/5 to 1/3 of that cost depending upon 
the number of stream tubes. The cost is roughly proportional to the chamber 
length analyzed (or the number of Az's), and to the number of stream tubes. Other 
influential variables are the numbers of spray group sizes, the number of 
corrector cycles utilized, the number of CRT plots and the iteration 
convergence tolerances . 

Thus, a complete DER analysis may vary in cos* om about 4 to 11 or more 

billing units. The example case given in Volume II cost 4.00 BU. 

Printout linecoui't depends upon the numbers of /Iz's, stream tubes, spray 
groups, iterations and print intervals. It is not unusual for one case 
print 10,000 lines. Concerning CRT's, rarely does one ask for more than 
10-12 graphs. 

The DER program requires two data set units other than the standard input/ 
output units. These units, which may be either magnetic tapes or disks, are 
used to transfer data between subroutines, both between and within major sub- 
program blocks. Usage is shown in Table 4. 


Table 4. Special Data Set Usage 


Datj 

e Set 

Usage 

No. 

Var. Name 

2 

M 

(1) Transfer Data from LISP to STC (STAPE) 

2 


(2) Save Data in STC for Plotting in STCRT 

3 


(1) Save and Retrieve Data for Iterating in STC 



(ITERS) 


3.5 PROGRAM OUTPUT 

Each of the major subprogram blocks has its own distinct output and each is 
discussed separately below. Print-out is the principal output and there are 
also some graphical and punched card outputs. 
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3>5tl LISP Output 

The output of the LISP computer program is provided both as the usual tabular 
printout, as computer^plotted CRT graphs and as a scratch tape record of data to 
be read and used by the STC computer program. 

A sample of LISP tabular output is presented in Volume II. First, there is a 
tabulation of all input data, vhich permits both a full documentation of the 
computer run conditions for later analysis and a convenient method to check 
input for errors if unusual results are calculated. The input data table is 
folloved by a one page table of specific spray distribution coefficient values 
used for each element specification. As in the example case, this table of 
coefficients rriay be preceded by a vaming message that some unlike doublet 
design or flow parameter falls outside the < .>rrelation range in the LISP library 
of coefficients along with values of pertinent parameters. None of the other 
coefficient subroutines provides such a message* 

Additional diagnostic date are printed out on the page following the distribu- 
tion coefficient values. These are: the sums of the elemental injection rates, 

SVEl and 2; the integrated sums of propellant flowrates calculated at the 
collection plane, TMCSl and 2; and, inverse collection efficiency factors, 
RC^NTl and 2, obtained by dividing SVEl by TMCSl , for example. Here, the 
1 and 2 refer to the fuel and oxidizer, respectively. Overall continuity is 
subsequently enforced by multiplying all mesh point flows and fluxes by BCjZlNTl 
or 2 , as appropriate. 

The input data is followed by a second table which cross references (by 
injector element) the calculated flow rates and drop sizes before evaporation 
to the read-in element coordinates 


The elesient reference table is followed by two tables referenced to the combustion 
zone mesh points. For the usual liquid/liquid injector, the third table lists 
the coordinates of the mesh points in the chamber slice at the collection plane 
z^. together with the weight flux (lb ^inch sec), the total collected mass 
(lb ^sec), the three mean droplet velocity components (ft/sec), and the mean 
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drop diameters (inch) of each propellant at the mesh points. The total 
collected mass at a mesh point is defined as the weight flux times the 
associated area at the mesh point. The values in this third table are based 
upon cold flow conditions, i.e., no vaporisation is assumed between elements 
and mesh point. The mesh points are listed in ascending order according to 
radial and angular coordinate. The last column of this table lists the sum of 
the collected mass of each propellant at all the mesh-points of constant radial 
coordinate, i.e., lists the radial distribution of the spray mass flux. 

If the injector uses gas/liquid elements, the definition of "Propellant 1" 
at the chamber mesh points is changed to mean the combustion gas resulting 
from the mixing of the injected gaseous fuel and the evaporated liquid 
oxidizer. Under these circumstances, the third table again lists the mesh 
point coordinates and local weight fluxes and collected mass of each propel- 
lant (Propellant 1 having its altered meaning), but the three spray velocity 
components and mean drop diameter are defined only for Propellant 2, liquid 
oxidizer. Listed in the place of the mean droplet diameter for Propellant 1 
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is the weight flux of evaporated oxidizer which arrives at the mesh point as 
part of the combustion gas. No listing of the collected propellant mass by 
radial coordinate Is made with gas liquid injectors. 

The fourth table again lists the coordinates of the mesh points in the chamber 
slice in the plane together with the reduced weight fluxes and droplet 
diameters of the collected spray after evaporation. A mass-weighted average 
evaporation of the original spray flux to each mesh point is also listed in 
this table. At the bottom of the fourth table are listed the Rupe mixing 
factor, the mixing limited c* efficiency, the overall percent vaporization 

of each propellant, and, finally, the average mixture ratio, terperaturc, 
velocity, and density of the (uniform) combustion gas at plane Zq. 

The fifth table, printed only for gas/liquid injection systems, again lists 
the mesh point coordinates, together with the associated liquid spray flux 
and drop diameter after evaporation, the fraction of the total oxidizer 
arriving at the mesh point which is in the vapor phase, the local gas mixture 
ratio, temperature and axial velocity, and the axial liquid spray velocity. 

A number of entries in this table represent repetitions of data in other 
tables which are grouped here with local gas variables for convenience in 
interpretation. 

The LISP graphical output is exemplified in Fig. 10 through lA. Fig. 10 
shows the mesh system for the chamber slice analyzed and the element origin 
locations for all injection elements considered to contribute flux to that 
slice. Fig. 11 is an example of the fuel and oxidizer mass flux profiles 
around the chamber slice at one fixed chamber radius. Figures 12 and 13 
are the contour plots of fuel and oxidizer mass flux for the entire chnmber 
cross-section. A similar plot for total mass flux is not shown. Figure 14 
is a contour plot of a modified fuel fraction function. The expression 
plotted is given at the top of the figure; it was chosen because it is 
bounded between zero and unity and has a value of 0.5 at the injection mixture 
ratio. 
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Figure 11. Fue^. and Oxidizer Mass Flux Profiles Computed 
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Figure 14. Contour Plot of Modified Fuel Fraction Computed by LISP 
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3,5«2 STC Output 

A sample case of STC computer program printout is included in Volume 11. Input 
data are vritten out imnediately as they are read in* This documents the data 
used for the particular case as veil as shoving vhether or not the data vere 
read^in properly* The input section should be examined for each case run to be 
sure that the intended input data vere actually used* 

Input data trasnferred from LISP are not printed out^ but a table of diagnostic 
data from subroutines STATE and SCBMBL is printed* Parameters appearing in tliat 
table are, vith flowrates in Ib/sec: 

Total fuel 9 oxidizer spray flovrates summed over all mesh 
points in STATE 

- Gaseous fuel, oxidizer flovrates 

- Total fuel, oxidizer flovrates transferred from LISP 

- The ratio TIFl/(SUMhl + WCF), etc. 

- The number of circular rings of mesh points from LISP 

- Total Flovrate assigned to the vail boundary layer strecun 
tube 

- The ratio 9fflL/(TIFl + TIF2) 

- The product PCTBL*(TIF1 + TIF2) 

- The indices of circular rings of mesh points in a given 
geometric zone 

- The cumulative total flovrate in a geometric zone and 
those set up prior to it 

The stream tube initialization data are tabulated and simultaneously punched 
out in cards* 

Based on the stream tube initialization data, a table is printed out of stream 
tube total flovrates and overall mixture ratios. This table is followed by 
values of the Rupe mixing efficiency factor, E^, and a mixing c'’^ efficiency for 
the stream tube flows* The latter represents an upper limit for multiple stream 
tube c*^ efficiency, since it corresponds to complete evaporation and burning of 
sprays within all stream tubes. 


SGMtfl, 2 

WGF, ^ 
TIFl, 2 
FP, 

K 

SffiL 

PCTT 

J^TT 

N1.N2 

SOM 
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Subroutine AYiUl sets up the array of chamber areas and writes out a table of 
chamber geometry information. Similar ly^ subroutine KI^IME computes and writes 
out tables of evaporation coefficients. 

Single stream tube analysis is preceded by writing out a one^page table of input 
total flows and averaged spray and gas parameters. During single stream tube 
analysis, data are written out as they are generated* At each z«-plane to be 
printed, complete gas and propellant spray group data are given. Additionally, 
the percentages of propellants evaporated and burned are listed and volume 
number mean propellant droplet diameters, computed* 

Calculated values of local flow Mach number and stagnation pressure (obtained 
from static pressure, Mach number and frozen ganma via the isentropic relation- 
ship) are printed in the right hand margin* In the printout of Volume II, 
these have been tom off* 

Two values each of flow area and contraction ratio may be given* Where the 
gas flow is subsonic, these should agree with each other precisely, whereupon 
the second set is not printed; if they do not there is some discrepancy between 
the read«*in AFR^(J,L) and AREAi(NT) arraya* At or near the nozsle throat 
plane, the two seta may disagree because the gas velocity has been set equal 
to sound speed and the local nozzle area adjusted to satisfy mass flow 
continuity* The contraction ratio calculated from continuity at the throat 
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plane is used as a multiplier (if it differs from unity by more than CRT0L) 
to adjust initial plane chamber pressure for a next iteration of single stream 
tube analysis. After the throat plane data are writte: out, a calculated c’*^ 
efficiency is listed. 

When the foregoing analysis has converged on its solution, the input value of 
nozzle radius ratio and calculated value of mean nozzle expansion coefficient, 
Y, are used by TRANS to generate transonic flow region isobars. The reduced 
coordinates and flow directions for each of 20 points along each isobar are 
written out, be^'^nning with the furthest downstream isobar and progressing 
upstream. Additionally, for the a*0 isobar, the absolute coordinates are 
written out for 40 points. Finally, a value is printed out for the nozzle 
discharge coefficient, 

TRANS also generates a CRT plot of the isobars* coordinates. Examples were 
shown earlier in Fig. 3. 

Multiple stream tube analysis follows the foregoing single stream tube and 
TRANS analyses. Stream tube input data are re-initialized and some additional 
data are written out to more completely define the initial-plane conditions. 
Initial-plane pressure is taken as the product of its value left in storage 
from the last preceding single stream tube iteration multiplied by the mixing 
c* efficiency and divided by TRANS *s nozzle discharge coefficient. 

At each prescribed z-plane for printing multiple stream tube results, complete 
definitive data for combustion gases and propellant sprays are written out. 
Local chamber area and contraction ratio are given; additionally, overall 
percentages of the propellants evaporated and burned are listed. Calculated 
stream tube stagnation pressures are printed in the right hand margin, which 
has been torn off in the printout of Volume II. 


At the throat position and intermittently downstream, diagnostic-tjrpe printouts 
containing data concerning dividing streamline intersections with the fifth, a«0, 
isobar are Inserted between the regular z-plane printouts. A summary table 
of these data is given near the end of the multiple stream tube printout. 

Finally, a long summary table is given of the stream tubes* outer radii at 
each z-plane. This is terminated with the minimum value of the sum of stream 
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tube areas and the ratio of that value to throat area ratio. The magnitude 
of this latter value’s deviation from unity determines whether or not all, 
or a portion*, of SIC’s multiple stream tube analysis will be repeated. If 
so, it is readily apparent in the printout. 

One (or more) computer-plotted graphs accompany the printout discussed above, 
showing the stream tubes’ outer radii along the entire chamber length, 
illustrated in Fig. 15. One of these graphs is plotted for each iteration 
through the multiple stream tube analysis. 

Following the last pass through the multiple stream tube analysis, a table 
is prin* d which lists the NAMELIST data punched by STC for subsequent use 
in running the improved TMC computer program. These are in the form of a TWC 
data deck, with all appropriate title cards, a problem definition card and 
NAMELIST data file cards. Only a part of the required data can be provided by 
STC, however, so comment cards are included to indicate the additional input 
data required before TDK can be run. STC provides values for the following 
narameters calculated at the tliroat plane: 


NZji^NES The number of stream tubes 

p(l) Stream- tube-area-weighted mean stagnation pressure 

XM(J) Mass fraction of the total gasified propellant flow 

carried by each stream tube 

^SKED(j)- The gasified propellant mixture ratio of each stream 

tube 

Finally, a table is printed out of additional data in anticipation of their 


being needed and used in a future revision of THC. Included are each stream 
tube’s cross-sectional area, ASUBS(j), static pressure, PSUBS(j) and stream gas 
velocity, VSUBS(j) at the position Z=ZPVSR where the TRANS pressure distribution 
is first invoked. Additional nozzle throat plane data are each stream tube’s 
ID’ “9 fraction (gas and spray) of the total injected propellant flowrate, XMl(j), 
overall mixture ratio, j^Fl(j), stagnation pressure, PNS(j), mean evaporation 

efficiency, ETAVAP(j), and mean fuel and oxidizer spray velocities, VBARF( j) and 
VBARi(j). 

^Subroutine ITER8 writes all pertinent data on scratch data unit 3 when the 
multiple stream tube analysis reaches the plane where absolute pressures 
are imposed on the nozzle. If only the subsequent portion of the analysis is 
to be iterated upon, ITER8 is recalled to read that record and reinitialize 
the analysis in that plane. 
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Figure 15, Example of Computer-Plotted Dividing 
Stream Lines from STC Program Block 
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